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ABSTRACT 

In this paper, we present and analyse optical photometry and spectra of the extremely 
luminous and slowly evolving Type la supernova (SN la) 2009dc, and offer evidence 
that it is a super-Chandrasckhar mass (SC) SN la and thus had a SC white dwarf (WD) 
progenitor. Optical spectra of SN 2007if, a similar object, arc also shown. SN 2009dc 
had one of the most slowly evolving light curves ever observed for a SN la, with a rise 
time of ~23 days and ATOi5(i?) = 0.72 mag. We calculate a lower limit to the peak 
bolometric luminosity of ~2.4 x 10*^ erg s~^, though the actual value is likely almost 
40% larger. Optical spectra of SN 2009dc and SN 2007if obtained near maximum 
brightness exhibit strong C II features (indicative of a significant amoimt of unburned 
material), and the post-maximum spectra are dominated by iron-group elements. All 
of our spectra of SN 2009dc and SN 2007if also show low expansion velocities. However, 
we see no strong evidence in SN 2009dc for a velocity "plateau" near maximum light 
like the one seen in SN 2007if (Scalzo ct al. 2010). The high luminosity and low 
expansion velocities of SN 2009dc lead us to derive a possible WD progenitor mass 
of more than 2 Mq and a ^^Ni mass of about 1.4-1.7 Mq. We propose that the host 
galaxy of SN 2009dc underwent a gravitational interaction with a neighboring galaxy 
in the relatively recent past. This may have led to a sudden burst of star formation 
which could have produced the SC WD progenitor of SN 2009dc and likely turned 
the neighboring galaxy into a "post-starburst galaxy." No published model seems to 
match the extreme values observed in SN 2009dc, but simulations do show that such 
massive progenitors can exist (likely as a result of the merger of two WDs) and can 
possibly explode as SC SNc la. 

Key words: supernovae: general - supernovae: individual (SN 2009dc, SN 2007if, 
SN 2003fg, SN 2006gz) 



1 INTRODUCTION 

Type la supernovae (SNe la) are differentiated from other 
types of SNe by the absence of hydrogen and the presence 
of broad absorption from Si II A6355 in tiicir optical spectra 
(for a review see Filippenko 1997). SNe la have been used 
to measure cosmological parameters to high precision (e.g., 
Kowalski et al. 2008; Hicken et al. 2009; Kessler et al. 2009; 
Amanullah et al. 2010), as well as to discover the accelerat- 
ing expansion of the Universe (Ricss ot al. 1998; Porlmuttcr 
et al. 1999). Broadly speaking, SNe la are the result of ther- 
monuclear explosions of C-O white dwarfs (WDs) resulting 
from either the accretion of matter from a nondegenerate 
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companion star (e.g., Whelan & Iben 1973) or the merger 
of two degenerate objects (e.g., Iben & Tutukov 1984; Web- 
bink 1984). However, the nature of the companion and the 
details of the explosion itself arc both still quite uncertain. 

The cosmological utility of SNc la comes from the fact 
that they are standardizable candles (i.e., their luminosities 
at peak can be calibrated). This naively seems reasonable 
since SNe la should all have the same amount of fuel and 
the same trigger point: they should all explode when a WD 
nearly reaches the Chandrasekhar mass of ~1.4 Mq. 

In 2003, SNLS-03D3bb (also known as SN 2003fg) was 
discovered (Howell ct al. 2006) and was shown to be a SN la 
that was overluminous by about a factor of 2, and had a 
slowly declining light curve, quite low expansion velocities, 
and unburned material present in near-maximum light spec- 
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tra. This last observation implies that a layer of carbon and 
oxygen from the progenitor existed on top of the burned sil- 
icon layer. In SNe la, carbon is usually extremely weak or 
completely absent (even at very early times) in optical and 
infrared (IR) spectra, although it has been (sometimes ten- 
tatively) identified in a few other cases (e.g.. Branch et al. 
2003; Marion et al. 2006; Thomas et al. 2007; Foley et al. 
2010). Howell et al. (2006) suggest that all of the oddities 
seen in SN 2003fg could be explained if its progenitor WD 
had a mass greater then the canonical upper limit for WDs 
of ~1.4 Mq, a so-called "super-Chandrasekhar mass" (SC) 
WD. 

Hickon et al. (2007) then presented data on SN 2006gz 
which shared some of the strange properties of SN 2003fg, 
and they concluded that SN 2006gz must have come from a 
WD merger leading to a SC SN la. Recently, Scalzo et al. 
(2010) published observations of SN 2007if which is yet an- 
other example of this emerging class of possible SC SNe la. 
They not only observed low expansion velocities, but they 
also saw a plateau in the expansion velocity near maximum- 
brightness which they interpret as evidence for the SN ejecta 
running into a shell of material. The accurately determined 
total mass of the SN 2007if system is well above the Chan- 
drasekhar mass. 

Recently, it was pointed out that there might be an- 
other member of this SC SN la class, SN 2009dc (Haru- 
tyunyan et al. 2009; Marion et al. 2009; Yamanaka et al. 
2009). SN 2009dc was discovered 15'.'8 west and 20'.'8 north 
of the nucleus of the SO galaxy UGC 10064 by Puckett et al. 
(2009) on 2009 Apr. 9.31 (UT dates are used throughout 
this paper), though in §3.1 wo will show a detection of the 
SN ~5 days earlier. It is located at aj2ooo = 15''51™12!12 
and i5j20oo ~ -|-25°42'28'.'0; the SN and its host are shown in 
Figure 1. No object was visible in our data at the position 
of the SN on 2009 Mar. 28 to a Umiting magnitude of ~19.5, 
so the actual explosion date was almost certainly between 
Mar. 28 and Apr. 4. 

Harutyunyan et al. (2009) obtained a spectrum of 
SN 2009dc one week after the announced discovery date 
and showed it to be a SN la before maximum light. They 
also noted that SN 2009dc spcctroscopically resembled the 
SC SN la candidate SN 2006gz (Hicken et al. 2007) at this 
time, but with a much lower expansion velocity as derived 
from the Si II A6355 absorption feature. Three days later, 
nearly simultaneous optical and IR spectra were obtained 
by Marion et al. (2009), covering a wavelength range of 0.36 
to 1.3 /xm. They again note some similarities to (as well as 
differences from) SN 2006gz, as well as similarities with the 
possible SC SN la 2003fg (Howell et al. 2006). 

Yamanaka et al. (2009) presented early-time optical 
and near-IR observations of SN 2009dc and showed that 
it did indeed share many of the properties seen in both 
SN 2003fg and SN 2006gz, suggesting that it too is a possi- 
ble SC SN la. Furthermore, Tanaka et al. (2010) published 
spectropolarimetry of SN 2009dc which indicated that the 
explosion was quite spherically symmetric. 

In this paper we present and analyse our own optical 
photometric and spectroscopic data for SN 2009dc (as well 
as spectra of SN 2007if) with the goal of more definitively 
answering the question of whether SN 2009dc and SN 2007if 
were truly SC SNe la. We show some of the most precise 
data on the rise time of a possible SC SN la ever published. 
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Figure 1. KAIT image of SN 2009dc and its host galaxy, 
UGC 10064. The field of view is 6.7'x6.7'. SN 2009dc is labeled 
along with the comparison stars used for differential photometry. 

The scale of the image is marked in the top right; north is up 
and east is to left. The SN is sufficiently far from its host galaxy 
that template subtraction was not performed when conducting 
photometry. 

as well as some of the latest-time photometric and spectral 

observations of a member of this class. In §2 we describe 
our observations and data reduction, and in §3 we discuss 
our analysis of the photometry and spectra of SN 2009dc 
(and its host galaxy) and SN 2007if. We attempt to robustly 
calculate physical parameters of SN 2009dc using a variety of 
methods and compare them to theoretical predictions in §4. 
Finally, in §5, we summarize our conclusions and ruminate 
about the future. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Photometry 

Observations of SN 2009dc began on 2009 Apr. 17, about 
one week before maximum B-band brightness, in BVRI fil- 
ters using the 0.76-m Katzman Automatic Imaging Tele- 
scope (KAIT; Filippenko et al. 2001) and the 1-m Nickel 
telescope, both at Lick Observatory. We continued to fol- 
low SN 2009dc for over 5 months until 2009 Sep. 26, when 
it reached the western limit of both telescopes in the early 
evening. We obtained late-time gVRI images of SN 2009dc 
using the dual-arm Low-Resolution Imaging Spectrometer 
(LRIS; Oko et al. 1995) with the 10-m Keck I telescope 
on 6 Feb. 2010 (281 days past maximum), BRI images us- 
ing the DEIMOS spectrograph (Faber et al. 2003) mounted 
on the 10-m Keck II telescope on 12 June 2010 (403 days 
past maximum), and F-band images again using LRIS on 
13 June 2010. 

Our optical photometry is complemented with data 
taken from the UltraViolet Optical Telescope (UVOT) on 
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the Swift Observatory in the U, B, V, UVWl, UVM2, and 
UVW2 filters. We downloaded 7 epochs of observations from 
the Swift archives. 

Our early-time optical images were reduced using 
a pipeline developed for KAIT and Nickel data (Gane- 
shalingam et al. 2010a). The images were bias subtracted 
and flatfielded at the telescope. Given the distance of 
SN 2009dc from the host galaxy, we did not find it necessary 
to attempt galaxy subtraction. We performed differential 
photometry using point-spread function (PSF) fitting pho- 
tometry on SN 2009dc and several comparison stars in the 
field with the DAOPHOT package in IRAF.^ The output 
instrumental magnitudes were calibrated to the standard 
Johnson BV and Cousins RI system using colour terms de- 
rived from many photometric nights. The comparison stars 
were calibrated against Landolt (1992) standard stars using 
13 photometric epochs to achieve a root-mean square (rms) 
of < 0.01 mag. The comparison stars arc identified in Fig- 
ure 1 with corresponding photometry in Table 1. To estimate 
the uncertainty in our photometry pipeline, we add artifi- 
cial stars with the same magnitude and PSF as the SN to 
each data image at positions of similar background bright- 
ness; the rms of doing this procedure 20 times is taken as 
the uncertainty. The photometric uncertainty is added in 
quadrature to the calibration error to produce our final un- 
certainty, adopting an error floor of 0.03 mag for B and / 
and 0.02 mag for V and R to account for slight systematic 
differences between the KAIT and Nickel photometry. Our 
final photometry of SN 2009dc is presented in Table 2. 

Late-time data obtained at the Keck I and Keck II tele- 
scopes were bias subtracted and flatfielded using standard 
imaging techniques. Differential photometry was performed 
using PSF fitting photometry on the SN and comparison 
stars that were not saturated, but also detected in our cal- 
ibration images obtained with the Nickel telescope. Cali- 
brations for the g band were obtained using the transfor- 
mations presented by Jester et al. (2005). In cases where 
all of the field stars from our Nickel calibration were satu- 
rated, calibrations for fainter stars in the field were obtained 
from the Sloan Digital Sky Survey (SDSS) and transformed 
into BVRI using transformations for stars from Jester et al. 

(2005) . Colour-term corrections were not applied. We in- 
clude a systematic error of 0.03 mag in all bands. Our final 
late-time photometry of SN 2009dc is presented in Table 3, 
which also includes an epoch of photometry obtained with 
the Nickel. 

We downloaded the Level-2 UVOT data from the Swift 
archive. Images taken during the same pointing were regis- 
tered and stacked to produce deeper images. We performed 
aperture photometry using the recipes prescribed by Li et al. 

(2006) for the optical data and Poole et al. (2008) for the UV 
data. We modified the ?7-band data to be in the Johnson- 
Cousins system using the colour corrections found in Li et al. 
(2006). In general, the UVOT B and V photometry is in 
good agreement with the photometry from our ground-based 



^ IRAF: The Image Reduction and Analysis Facility is dis- 
tributed by the National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in As- 
tronomy (AURA) under cooperative agreement with the National 
Science Foundation (NSF). 
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Table 5. Journal of Spectroscopic Observations of SN 2009dc 



UT Date Age^ Range (A) Airmass"^ Exp (s) 



2009 Apr. 18.5 


-7 


3500- 


-9900 


1.06 


1500 


2009 May 31.3 


35 


3500- 


-10200 


1.03 


1500 


2009 Jun. 17.5'= 


52 


3400- 


-10200 


1.56 


250/200=^ 


2009 Jun. 29.3 


64 


3500- 


-10200 


1.08 


2100 


2009 Jul. 16.3 


80 


3500- 


-10200 


1.29 


2100 


2009 Jul. 23.2 


87 


3500- 


-10200 


1.10 


2100 


2009 Jul. 28.3 


92 


3500- 


-10200 


1.33 


2100 


2009 Aug. 14.2 


109 


3500- 


-10000 


1.31 


2400 


2010 Feb. 6.6" 


281 


3500- 


-10200 


1.19 


3 X (630/600)'* 



^Rest-frame days relative to the date of B-band brightness max- 
imum, 2009 Apr. 25.4 (sec §3.1). 
''Airmass at midpoint of exposure. 

'^Thcsc observations used LRIS (Okc et al. 1995) on the 10-m 
Keck I telescope. The others used the Kast spectrograph on the 
Lick 3-m Shane telescope (Miller &c Stone 1993). 
''The blue side was exposed longer than the red side due to the 
relatively long readout time of the red-side CCD in LRIS. 

telescopes to within ±0.05 mag. The UVOT S-band data are 

systemically brighter, which could possibly be attributed to 
galaxy light falling within our aperture. Our results are given 
in Table 4. 

We did not follow the photometric behaviour of 
SN 2007if, another SC SN la candidate, but we do present 

spectroscopic observations (see §2.2). Throughout the rest of 
the paper we adopt 2007 Sep. 5.4 as the date of B-band max- 
imum brightness and 0.07416 as the redshift (z) of SN 2007if, 
both taken from Scalzo et al. (2010). 



2.2 Spectroscopy 

Beginning about a week before maximum brightness, optical 
spectra of SN 2009dc were obtained mainly using the dual- 
arm Kast spectrograph (Miller & Stone 1993) on the Lick 
3-m Shane telescope. Two spectra were also obtained us- 
ing LRIS on Keck I. Our last spectral observation occurred 
281 days after B-band maximum. 

The Kast spectra all used a 2" wide slit, a 600/4310 
grism on the blue side, and a 300/7500 grating on the red 
side, yielding full-width at half-maximum (FWHM) resolu- 
tions of ~4 A and ~10 A, respectively. The LRIS spectrum 
was obtained with a 1" slit, a 600/4000 grism on the blue 
side, and a 400/8500 grating on the red side, resulting in 
FWHM resolutions of ~4 A and ~6 A, respectively. All ob- 
servations were aligned along the parallactic angle to reduce 
differential light losses (Filippenko 1982). Table 5 summa- 
rizes the spectral data of SN 2009dc presented in this paper. 

We obtained three spectra of SN 2007if using LRIS. For 
the first two LRIS spectra, we used a 400/3400 grism on the 
blue side and a 400/8500 grating on the red side, resulting 
in FWHM resolutions of ~ 6 A on both sides, while the 
final spectrum of SN 2007if was obtained with a 600/4000 
grism on the blue side, giving a resolution of ~4 A. In all 
cases, the long, l"-widc slit was aligned along the parallactic 
angle to reduce differential light losses. Table 6 summarizes 
the spectral data on SN 2007if presented here. We also note 
that our last spectrum of SN 2007ff (from 2007 Dec. 13) 
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Table 1. Comparison Stars for the SN 2009dc Field 



Star 


Q^J200(J 


,12000 


B (mag) 


V (mag) 


R (mag) 


/ (mag) 




SN 




-l^9^-zL9-95^ n 












1 


15:51:00.40 


+25:44:00.1 


18.769 (010) 


17.418 (006) 


16.503 (005) 


15.744 (Oil) 


9 


2 


15:50:59.16 


+25:43:04.7 


17.501 (008) 


16.957 (006) 


16.608 (004) 


16.236 (008) 


12 


3 


15:51:00.72 


+25:43:02.5 


18.126 (009) 


16.816 (004) 


15.940 (004) 


15.228 (008) 


10 


4 


15:51:16.96 


+25:42:40.6 


18.671 (006) 


17.805 (005) 


17.249 (006) 


16.734 (Oil) 


11 


5 


15:51:11.83 


+25:41:48.3 


16.661 (007) 


15.872 (004) 


15.406 (003) 


14.984 (009) 


13 


6 


15:51:13.64 


+25:41:20.5 


15.383 (007) 


14.471 (004) 


13.935 (003) 


13.460 (009) 


13 


7 


15:51:01.38 


+25:41:11.2 


15.892 (007) 


15.200 (005) 


14.808 (003) 


14.420 (010) 


13 


8 


15:51:17.44 


+25:39:59.6 


17.156 (006) 


16.534 (004) 


16.174 (004) 


15.790 (008) 


13 


1(7 uncertainties are in 


parentheses, 


in units of 0.001 


mag. 









Table 6. Journal of Spectroscopic Observations of SN 2007if 



UT Date 


Age'' 


Range (A) 


Airmass'' 


Exp. (s) 


2007 Oct. 15.4 


37 


3300-9200 


1.03 


1200 


2007 Nov. 12.2 


63 


3300-9150 


1.49 


1200 


2007 Dec. 13.2^^ 


92 


3300-9150 


1.01 


1200 



^Rest-frame days relative to the date of B-band maximum bright- 
ness, 2007 Sep. 5.4 (Scalzo et al. 2010). 
''Airmass at midpoint of exposure. 

'^This observation used a slightly higher resolution grism for the 
blue side. 



was taken under nonideal observing conditions (clouds were 
present and the atmospheric seeing was poor), and thus its 
spectrophotometric accuracy is not as good as that of the 
other observations. 

All spectra were reduced using standard techniques 
(e.g., Foley ot al. 2003). Routine CCD processing and spec- 
trum extraction were completed with IRAF, and the data 
were extracted with the optimal algorithm of Horne (1986). 
We obtained the wavelength scale from low-order polyno- 
mial fits to calibration-lamp spectra. Small wavelength shifts 
were then applied to the data after cross-correlating a tem- 
plate sky to the night-sky lines that wore extracted with 
the SN. Using our own IDL routines, we fit spectrophoto- 
metric standard-star spectra to the data in order to flux 
calibrate our spectra and to remove telluric lines (Wade & 
Horne 1988; Mathcson et al. 2000). Information regarding 
both our photometric and spectroscopic data (such as ob- 
serving conditions, instrument, reducer, etc.) was obtained 
from our SN database (SNDB). The SNDB uses the popu- 
lar open-source software stack known as LAMP: the Linux 
operating system, the Apache webserver, the MySQL rela- 
tional database management system, and the PHP server- 
side scripting language (see Silverman et al. 2010, for further 
details). 



3 ANALYSIS 
3.1 Light Curves 

We present our final optical BVRI light curves of SN 2009dc 
in Figure 2. We include for comparison a "normal" Type la 
SN 2005cf (Wang et al. 2009b), another SC SN la candidate 
SN 2006gz (Hicken et al. 2007), the "standard overlumi- 



nous" Type la SN 1991T (Lira et al. 1998), and the peculiar 

SN 2002cx-hke SN 2005hk (Phillips et al. 2007). 

The light curves of SN 2009dc have many features 
which are noticeably distinct from most other SNe la. The 
light curves are much broader than those of spectroscopi- 
cally normal SNe la such as SN 2005cf. Even in comparison 
to the ovcrluminous SN 1991T, SN 2009dc evolves signifi- 
cantly more slowly. The light curve of SN 2006gz presented 
by Hicken et al. (2007) provides an good match, although 
SN 2009dc appears to have a slower rise time and a slower 
decline. 

The absence of a prominent second maximum in the 
R and / bands is particularly striking. The secondary 
maximum in I has been attributed to the propagation of 
an ionization front through iron-group elements (IGEs) as 
they transition from being doubly ionized to singly ionized 
(Kasen 2006). Observationally, the strength of the secondary 
maximum has been found to correlate with peak luminos- 
ity, powered by Ni decay, with luminous SNe exhibiting more 
prominent secondary peaks (Hamuy et al. 1996; Nobili et al. 
2005). However, this feature is weak in SN 2009dc, even 
though it is more luminous than a typical SN la. Kasen 
(2006) found that significant mixing of ®^Ni in the compo- 
sition of a SN la can lead to a blending of the two J-band 
peaks, possibly explaining why a secondary peak is weak in 
SN 2009dc. 

After correcting for Milky Way (MW) extinction (from 
the dust maps of Schlegel et al. 1998), we fit a polynomial 
to the B-band light curve to find that SN 2009dc peaked 
on JD 2,454,946.93 ± 0.2 (on 2009 Apr. 25.4) at B^ax = 
15.05 ± 0.04 mag. These values are consistent with those 
derived by Yamanaka et al. (2009) . 

We also measure Ami6(B), the decline in flux from 
maximum light to 15 days past maximum in the B band. 
Correcting for time dilation, we find Ami5(S) = 0.72 ± 
0.03 mag, making SN 2009dc one of the most slowly evolv- 
ing SNe ever discovered and comparable to other SC SN la 
candidates in the literature. Interestingly, Yamanaka et al. 
(2009) measure Ami6(B) = 0.65 ± 0.03 mag. Comparing 
the two photometric datasets, we find discrepancies between 
measurements of local standard stars of ~0.15 mag. With 
over 10 calibrations of the field on photometric nights, we 
are confident in our measurements of the field standards 
presented in Table 1. While these differences are troubling, 
they do not alter the conclusion that SN 2009dc is a slowly 
evolving SN with an extreme value of Ami5(B). 

The host galaxy of SN 2009dc is part of the Lick Obser- 
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Table 2. Early-time BVRI photometry of SN 2009dc 



JD B (mag) V (mag) _R (mag) / (mag) Telescope 



z454yz5 


95 














18 


824 


(177) 








KAIT 


245490S 


98 


15 


670 


(030) 


15 


638 


(020) 


15 


596 


(020) 


15 


718 


(030) 


KAIT 


24549^9 


95 


15 


549 


(030) 


15 


541 


(020) 


15 


500 


(020) 


15 


577 


(030) 


Nickel 


CI A C Art Af\ 

2454940 


92 


15 


474 


(030) 


15 


472 


(030) 


15 


430 


(030) 


15 


524 


(030) 


Nickel 


z4o4y4U 


93 


15 


oil 


[030) 


15 


495 


[020) 


15 


475 


[020) 


15 


608 


[030) 


TaT a TT" 
r\ Al i 


2454942 


95 


15 


455 


(030) 


15 


409 


(020) 


15 


376 


(020) 


15 


529 


(030) 


KAIT 


24o494d 


97 


15 


383 


(030) 


15 


327 


(020) 


15 


272 


(020) 


15 


450 


(030) 


i\ Al 1 


2454948 


87 


15 


369 


(030) 


15 


350 


(020) 


15 


288 


(020) 


15 


454 


(030) 


KAIT 


24o49ol 


92 


15 


449 


(030) 


15 


386 


(020) 


15 


311 


(020) 


15 


449 


(030) 


i\ Al 1 


2454957 


97 


15 


770 


(030) 


15 


535 


(020) 


15 


423 


(020) 


15 


537 


(030) 


KAIT 


2454909 


85 


15 


930 


[030) 


15 


610 


[020) 


15 


AitA 

4o4 


[020) 


15 


554 


[030) 


TV" A TT* 

K Al i 


2454950 


92 


16 


004 


(030) 


15 


591 


(020) 


15 


476 


(020) 


15 


593 


(030) 


KAIT 


24o49D2 


93 


16 


137 


[036) 


15 


656 


[024) 


15 


543 


[026) 


15 


633 


[084) 


jy A Trp 
K.Ai i 


CI A c Ana A 

24549b4 


90 


16 


322 


(030) 


15 


762 


(020) 


15 


584 


(020) 


15 


568 


(030) 


KAIT 


24549DO 


93 


16 


528 


(030) 


15 


846 


(020) 


15 


642 


(020) 


15 


562 


(030) 


KAIT 


24549DD 


95 


16 


478 


(030) 


15 


839 


(020) 


15 


611 


(020) 


15 


481 


(030) 


Nickel 


24549DO 


89 


16 


677 


(030) 


15 


926 


(020) 


15 


675 


(020) 


15 


565 


(030) 


KAIT 


245497U 


88 








16 


029 


(086) 


15 


696 


(053) 








KAIT 


2454971 


92 


16 


880 


(030) 


16 


054 


(020) 


15 


701 


(024) 


15 


500 


(030) 


Nickel 


c\ A^ A nTO 

2454972 


89 


16 


972 


(030) 


16 


106 


(020) 


15 


748 


(020) 


15 


579 


(030) 


KAIT 


2454974 


88 


17 


154 


(030) 


16 


202 


(020) 


15 


799 


(020) 


15 


557 


(030) 


KAIT 


2454975 


87 


17 


188 


(030) 


16 


239 


(020) 


15 


797 


(020) 


15 


526 


(030) 


Nickel 


n A i^A n'7'7 
2454977 


9U 


17 


384 


[030) 


16 


334 


[020) 


15 


OTA 

o74 


[020) 


15 


659 


[030) 


rvAi 1 


2454981 


89 


17 


650 


(047) 


16 


488 


(020) 


16 


009 


(020) 


15 


707 


(030) 


KAIT 


2454989 


83 


18 


045 


[078) 


16 


893 


[027) 
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Table 3. Late-time photometry of SN 2009dc 



JD 


Phase^ 




Filter 




Mag 






9^0 
zou 






uuu 


Z 1 .ouo 






250 




Y 


360 


21.016 


0.320 


2455233.10 


281 


Keck/LRIS 


9 


180 


21.894 


0.050 




281 


Keck/LRIS 


V 


180 


21.988 


0.042 




281 


Keck/LRIS 


R 


60 


22.600 


0.084 




281 


Keck/LRIS 


I 


120 


21.483 


0.060 


2455359.05 


403 


Kcck/DEIMOS 


B 


360 


25.010 


0.120 




403 


Kcck/DEIMOS 


R 


450 


24.987 


0.143 




403 


Keck/DEIMOS 


I 


450 


23.746 


0.185 


2455360.10 


404 


Keck/LRIS 


V 


540 


24.834 


0.152 



"■Rest-frame days relative to the date of B-band maximum brightness. 



Table 4. UVOT Photometry of SN 2009dc 



JD 


Filter 


Mag 


(7 


JD 


Filter 


Mag 




2454946.57 


UVW2 


17.338 


0.044 


2454946.57 


U 


14.720 


0.017 


2454951.85 


UVW2 


17.756 


0.079 


2454952.21 


U 


15.124 


0.010 


2454952.22 


UVW2 


17.987 


0.073 


2454955.72 


U 


15.506 


0.022 


2454955.72 


UVW2 


18.246 


0.053 


2454960.61 


u 


16.016 


0.027 


2454960.61 


UVW2 


18.741 


0.066 


2454980.35 


u 


17.890 


0.066 


2454980.35 


UVW2 


20.220 


0.146 


2454984.44 


u 


18.070 


0.088 


2454984.44 


UVW2 


20.314 


0.177 


2454946.57 


B 


15.324 


0.017 


2454946.58 


UVM2 


17.261 


0.027 


2454955.72 


B 


15.630 


0.017 


2454955.73 


UVM2 


18.346 


0.058 


2454960.61 


B 


15.935 


0.019 


2454960.61 


UVM2 


18.884 


0.075 


2454980.35 


B 


17.542 


0.042 


2454980.35 


UVM2 


19.986 


0.135 


2454984.44 


B 


17.739 


0.055 


2454984.45 


UVM2 


19.586 


0.138 


2454946.58 


V 


15.286 


0.028 


2454946.53 


UVWl 


16.041 


0.032 


2454955.73 


V 


15.489 


0.027 


2454951.71 


UVWl 


16.537 


0.035 


2454960.61 


V 


15.587 


0.028 


2454955.65 


UVWl 


16.965 


0.035 


2454980.35 


V 


16.550 


0.044 


2454960.27 


UVWl 


17.404 


0.042 


2454984.45 


V 


16.714 


0.058 


2454980.20 


UVWl 


18.644 


0.077 










2454984.17 


UVWl 


18.812 


0.099 











vatory Supernova Search (LOSS; Li et al. 2000; Filippenko 
et al. 2001), allowing us to put strict constraints on the rise 
time of SN 2009dc. Our first detection of SN 2009dc is from 
an unfiltered image on 2009 Apr. 4, when the SN is seen at 
R = 18.82±0.18 mag, indicating a rise time > 21 days. This 
first detection is included in Figure 2 as part of our i?-band 
data (Li et al. 2003, show that KAIT unfiltered data approx- 
imate the R band). In an image taken on 2009 Mar. 28, there 
is no detection of SN 2009dc down to R ^ 19.3 mag. With 
these constraints wc conservatively estimate a rise time of 
23 ± 2 days. Haydcn et al. (2010) recently found an average 
SN la rise time of 17.38±0.17 days using 105 SDSS-II SNe la, 
with slowly declining SNe tending to have shorter rise times. 
The rise time of SN 2009dc is significantly longer than the 
average SN la in their sample, despite being a slowly declin- 
ing SN, and therefore does not follow the trend found for 
normal SNe la in their sample. Riess et al. (1999), on the 
other hand, found an average rise time of 19.5 ± 0.2 days 
for a "typical" SN la (Ami5(B) = 1.1 mag), with slowly 
declining SNe having slower rise times. Their sample of 10 
objects lacks SNe la with Am(15)s < 0.95 mag, making an 
extrapolation to Ami5(S) = 0.72 mag to compare with the 
rise time of SN 2009dc unreliable. 

The cosmological application of SNe la as precise dis- 
tance indicators relies on being able to standardize their 
luminosity. Phillips (1993) showed that Ami6(B) is well cor- 



related with luminosity at peak brightness for most SNe la, 
the so-called "Phillips relation". In Figure 3, we plot abso- 
lute V magnitude as a function of Ami5(i?) for 71 SNe la 
with zvirgo infaii > 0.01 from the LOSS photometry database 
(Ganeshalingam et al. 2010b) and SN 2009dc to determine 
whether SN 2009dc follows the Phillips relation. Host-galaxy 
extinction for the sample of 71 SNe la is derived using 
MLCS2k2.v006 with the g/osz prior, which assumes that the 
late-time {B — V)t^+3z a colour is indicative of the host- 
galaxy extinction (Jlia et al. 2007). For SN 2009dc, the 
late-time colour evolution differs significantly from that of 
normal SNe la; thus, we instead adopt E{B — F)host = 
0.1 mag, assuming Ry = 3.1 for the host galaxy (see §3.5 
for details). To place the SNe on an absolute scale, we 
use the luminosity distance assuming a ACDM cosmology 
with Hq = 70 km s"^Mpc"^ The best-fit line to SNe with 
0.6 < Ami5(B) < 1.7 mag, excluding SN 2009dc, is plotted 
in Figure 3, with the area that falls within la of the relation 
shaded in grey. SN 2009dc is clearly an overluminous outlier 
in comparison to SNe la having similar values of Ami5{B). 
Events similar to SN 2009dc cannot be standardized using 
current distance-fitting techniques which rely on parameter- 
izations of light-curve shape (Howell et al. 2006). 

We measure the late-time decay of the light curves using 
a linear least-squares fit to data in the range ~50-150 days 
after maximum. We find a decline of 1.47 ± 0.04 mag per 
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Figure 2. Optical and ultraviolet (UV) light curves of SN 2009dc 
from KAIT and the 1-m Lick Nickel telescope. For comparison 
we also plot the optical light curves of the overluminous Type la 
SN 1991T, the normal SN 2005cf, the peculiar SN 2002cx-like 
SN 2005hk, and another SC SN la candidate SN 2006gz. The 
light curves of each SN are shifted so that t = days corresponds 
to the time of Bmax and reach the same peak magnitude for a 
given filter. The evolution of SN 2009dc is atypical compared to 
other SNe la. SN 2009dc evolves much more slowly and exhibits 
only a muted secondary maximum in the / band. We include 
in our _R-band light curve an early detection of SN 2009dc {red 
open circle) in an unfiltered image from the Lick Observatory SN 
Search (LOSS), about 3 weeks before maximum brightness. The 
data sources for the other SNe are cited in the text. 



100 days in B, 1.91 ± 0.02 mag per 100 days in V, 2.78 ± 
0.03 mag per 100 days in R, and 2.87±0.04 mag per 100 days 
in I. Leibundgut (2000) find typical decay rates for SNe la of 
1.4 mag per 100 days in B, 2.8 mag per 100 days in V , and 
4.2 mag per 100 days in I. SN 2009dc shows significantly 
slower decline rates in V and I. 

Late-time photometry of SN 2009dc (Table 3) taken 
with the Nickel telescope shows only marginal detections 
in B and V and upper limits in R and 7, while our Keck 
images, 281 days past maximum, give clear detections in 
gVRI. Figure 4 shows the late-time behaviour of SN 2009dc 
in comparison to late-time photometry of SN 2003du, a typ- 
ical SN la, taken from Stanishev et al. (2007), and the linear 
decay rates found using data ~50-150 days after maximum. 
A constant has been added to each SN 2003du light curve 
in order to match the peak magnitude of SN 2009dc. Com- 
pared to the extrapolations based on the measured linear 




1.0 \2 1.4 1.6 

A »ii5(B) [mag] 



2.0 



Figure 3. My as a function of Ami5(_B) for 71 SNe la 
with zvirgo infall > 0-01 from the LOSS photometry archive 
{filled circles, Ganeshalingam et al. 2010b), compared with 
SN 2009dc {filled star). Host-galaxy extinction is derived using 
MLCS2k2.v006 (Jha et al. 2007), except in the case of SN 2009dc 
where we adopt E{B — y)host =0.1 mag with Ry = 3.1 (see 
§3.5 for details). The luminosity distance was calculated using 
the ACDM concordant cosmology with = 70 km s~^Mpc~^, 
Um = 0.27, and Ua = 0.73 (Spergel et al. 2007). The best-fit line 
to the data, excluding SN 2009dc, in the range 0.6 < Ami5{B) < 
1.7 mag is plotted as a dashed line, with points that fall within 
Icr shaded in grey. SN 2009dc is clearly an outlier that does not 
follow the luminosity-width relation. 



decay rates, at 281 days past maximum SN 2009dc is fainter 
by ~ 0.5 mag in B and V, brighter by ~0.5 mag in R, 
and ~1.5 mag brighter in /. In comparison to SN 2003du, 
SN 2009dc is within ~0.15 mag of the interpolated values 
in BVR and ~0.5 mag brighter in I. We caution that our 
interpolated values for SN 2003du suffer from a rather large 
gap in the light curve between 225 days and 366 days af- 
ter maximum. While the late-time behaviour of SN 2009dc 
does not match the linear decay of SN 2003du exactly in 
all bands, our detections indicate that it has not undergone 
an unexpected drop in luminosity, and it is consistent with 
the late-time light curve being powered by ^^Co decay (see 
§4.2). 

The data taken ~400 days past maximum show only 
marginal detections in all bands, indicating that the SN is 
still active despite the points lying below interpolated values 
from SN 2003du. The detections in B and V are clearly 
below what we would expect from SN 2003du, while R and 
I are not too far below expectations. We address possible 
reasons for the drop in flux in §4.2.2. 



3.2 Colour Evolution 

The colour evolution in B — V, V — R, and R — I for 
SN 2009dc in comparison to SNe 2006gz, 2005cf, 1991T, 
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Figure 4. The late-time decay of SN 2009dc from 50 to about 
400 days past maximum light compared to that of SN 2003du 
{solid line), a typical SN la. The light curves are shifted such 
that t = days corresponds to the time of -Bmax, and a constant 
has been added to each SN 2003du light curve to match the peak 
magnitude of SN 2009dc. Filled circles are KAIT and Nickel data, 
whereas open squares are LRIS and DEIMOS data. Upper limits 
are marked with arrows. The LRIS g-band point is plotted in the 
B-band light curve panel. We fit the decay in BVRI, plotted as 
a dashed line, using a linear least-squares fit to our well-sampled 
data between 50—150 days past maximum, before the SN became 
projected too close to the Sun. Late-time photometry obtained 
with LRIS at 281 days past maximum indicates that the flux in 
B and V is below what is expected from linear extrapolations, 
but mostly agrees with expectations from SN 2003du. However, 
the flux in R and / is larger than what is expected compared to 
the extrapolations. Our _R-band detection falls on the comparison 
light curve, while the /-band detection is brighter than what is 
expected. SN 2009dc is detected in all bands at '~400 days past 
maximum, indicating that the SN is still active despite the points 
lying below interpolated values from SN 2003du. The data for 
SN 2003du were taken from Stanishev et al. (2007). 



and 2005hk is displayed in Figure 5. All SNe have been 
corrected for MW extinction using reddening derived from 
the dust maps of Schlegel et al. (1998). We have also cor- 
rected for host-galaxy extinction, assuming Rv = 3.1, using 
the following reported values of E{B — V)host: 0.13 mag for 
SN 1991T (Lira et al. 1998), 0.10 mag for SN 2005cf (Wang 
et al. 2009b), 0.09 mag for SN 2005hk (Phillips et al. 2007), 
0.18 mag for SN 2006gz (Hicken et al. 2007), and 0.10 mag 
for SN 2009dc. 

The colour curves indicate that SN 2009dc was a partic- 
ularly blue SN la even compared to the prototypical over- 
luminous SN 1991T. In B-V, SN 2009dc remains bluer 
than SNe 1991T, 2005cf, and 2005hk until ~50 days past 
maximum light. The colour curves for SN 2009dc are most 
similar to those of SN 2006gz, another SC SN la candidate. 
At i > 50 days, the B - V colour of SN 2009dc becomes 
redder than that of SN 1991T and SN 2005cf. 
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Figure 5. B-V (top), V-R (middle), and R-I (bottom) 
colour curves of SN 2009dc. Plotted for comparison are the colour 
curves of SNe 1991T, 2005cf, 2005hk, and 2006gz. All curves have 
been corrected for MW reddening and host-galaxy extinction. 
We have not included errors in derived host-galaxy extinction for 
SN 2009dc which will systematically shift curves in one direction. 
At early times, SN 2009dc is bluer than SNe 1991T, 2005cf, and 
2005hk. The evolution most clearly resembles that of SN 2006gz. 
We plot the Lira law as a solid line in the range 35 < t < 85 days. 
SN 2009dc shows a slower red to blue colour evolution compared 
to the Lira-Phillips relation, especially after t = 70 days. The 
data sources are cited in the text. 



Lira (1996) showed that SNe with low host-galaxy ex- 
tinction had similar B — V colour evolution between t — -f 30 
to -f 90 days independent of light-curve shape (the "Lira- 
Phillips relation"). Hicken et al. (2007) used this relationship 
to derive the host-galaxy extinction for SN 2006gz. However, 
a comparison between the slope of the Lira-Phillips rela- 
tion and the B — V colour evolution of SN 2009dc shows 
disagreement. Adopting the relation derived using 6 low- 
reddening SNe la by Phillips et al. (1999) and fitting for 
a constant E{B — y)host between t — +35 and -f 90 days, 
we find E{B - V)host = 0.26 ± 0.04 mag, with = 41 
for 15 degrees of freedom. If we restrict our fit to be- 
tween t = -1-35 and -f70 days, we find a better fit with 
E{B - V")host = 0.19 ± 0.04 mag; = 5 for 7 degrees of 
freedom. We include the fit to the Lira-Phillips relation in 
Figure 5. 
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3.3 Pre-Maximum Spectrum 

On 2009 Apr. 16.22, ~ 9 days before B-band maximum, 
Harutyunyan ct al. (2009) obtained a spectrum of SN 2009dc 
which showed "prominent C II lines and a blue continuum" , 
and they noted that the spectrum resembled pre-maximum 
spectra of SN 2006gz (Hickon ot al. 2007), but with a much 
lower expansion velocity. We obtained a spectrum ~2 days 
later which confirmed the spectral peculiarities noted by 
Harutyunyan et al. (2009). 

Our pre-maximum spectrum is shown in Figure 6, 
where wc compare SN 2009dc to other possible SC SNe la: 
SN 2006gz (Hicken et al. 2007), SN 2003fg (Howell et al. 
2006), and SN 2004gu (Contreras et al. 2010, though the 
displayed spectrum is from our own database). In addi- 
tion, we show for comparison spectra from similar epochs 
of the "standard overluminous" Type la SN 1991T (Filip- 
penko et al. 1992), the SN 2002cx-like peculiar SN 2005hk 
(Chornock et al. 2006), and the "standard normal" Type la 
SN 2005cf (Wang ct al. 2009b). 

SN 2009dc has strong Si II and S II lines, and appears 
to also contain O I, all of which are usually found in near- 
maximum spectra of SNe la. However, there are also two ap- 
parent C II lines (labeled in Fig. 6 and discussed in §3.3.2) 
seen in SN 2009dc. The existence (and strength) of these 
features is one reason that SN 2009dc was immediately clas- 
sified as a possible SC SN la (Harutyunyan et al. 2009), and 
one way used to distinguish it from other subtypes of SNe la. 

Spectroscopically, SN 2009dc does not seem to resemble 
any of SNe 1991T, 2005hk, or 2005cf. Although SN 2006gz 
has weak C II features at 14 days before maximum, by 8 days 
before maximum they have almost completely disappeared 
(whereas SN 2009dc has strong C II at 7 days before max- 
imum). In addition, it is apparent from Figure 6 that the 
expansion velocity of SN 2006gz is much larger than that 
of SN 2009dc (based on the Si II A6355 feature, for exam- 
ple; see §3.3.1 for further details). Finally, SN 2003fg does 
seem to resemble SN 2009dc, even though the spectrum of 
SN 2003fg has low signal-to-noise ratio (S/N) and the C II 
features are not nearly as prominent in SN 2003fg as they are 
in SN 2009dc. It is possible that this is due to the fact that 
the spectrum of SN 2003fg was taken 9 days later (relative 
to maximum light) than our spectrum of SN 2009dc. 



3.3.1 Sill 

The tell-tale spectroscopic signature of a SN la is broad 

absorption from Si II A6355 (e.g., Filippcnko 1997), and 
SN 2009dc is no different. One thing that is interesting about 
the Si II feature in SN 2009dc is its extremely low velocity. 
The minimum of the Si II A6355 absorption feature was 
found to be blueshifted by 8700 km s"'^ about 9 days be- 
fore maximum brightness (Harutyunyan et al. 2009), and 
we measure the feature to bo blueshifted by ~8600 km s~^ 
in our spectrum taken 7 days before maximum. Yamanaka 
et al. (2009) report it to be blueshifted by 8000 km s"^ 
in their spectrum 3 days before maximum, decreasing to 
6000 km s^^ by 25 days after maximum. 

These velocities are much lower than the average photo- 
spheric velocity near maximum for SNe la of ~11,000 km 
(Wang et al. 2009a). In fact, according to Wang et al. 



(2009a), the Si II velocities seen in SN 2009dc are approxi- 
mately 8a below the average.^ 

Interestingly, these velocities are similar to those of 
two other possible SC SNe la, SN 2003fg and SN 2007if 
(Howell et al. 2006; Scalzo et al. 2010, respectively). How- 
ever, two other objects that were claimed to be possi- 
ble SC SNe la show much more normal Si II velocities. 
Hicken et al. (2007) report expansion velocities near max- 
imum light of about 11,000-13,000 km s"^ for SN 2006gz. 
SN 2004gu, which was compared to SN 2006gz by Contr- 
eras et al. (2010), also shows a normal photospheric velocity 
of ~11,500 km at 5 days before maximum (this work) 
slowing to ~11,000 km s^^ at 3 days past maximum.^ 

Figure 4 of Yamanaka et al. (2009) and Figure 4 of 
Scalzo et al. (2010) both present the velocity evolution of 
the Si II A6355 feature for various possible SC SNe la and 
other comparison SNe. Our new early-time data point (i.e., 
~8600 km s~^ at f « —7 days) is quite consistent with the 
rest of the published SN 2009dc data. Moreover, we detect 
Si II A6355 in our next spectrum taken 35 days past max- 
imum at a blueshifted velocity of about 5500 km s~^ (see 
§3.4.1 for further details). This yields a velocity gradient of 
~74 km s^"*^ day""*^, which is over twice as large as that found 
for SN 2007if (34 km s"^ day"^ Scalzo ct al. 2010). Further- 
more, combining our measurements with all previously pub- 
lished values of the expansion velocity of SN 2009dc (Haru- 
tyunyan et al. 2009; Yamanaka et al. 2009; Tanaka et al. 
2010), wc sec no strong evidence for a velocity "plateau" 
near majcimum light like the one seen in SN 2007if (Scalzo 
et al. 2010). 

We measure the equivalent width (EW) of the Si II 
A6355 feature to be ~ 40 A in our spectrum of SN 2009dc 
obtained 7 days before maximum. This is .similar to the EW 
of the same line in SN 2006gz at similar epochs (Hicken 
et al. 2007), and both are well below the average Si II EW 
for normal SNe la (though they are comparable to the over- 
luminous SN 1991T-like SNe la, e.g., Hachinger et al. 2006; 
Wang et al. 2009a). 

3.3.2 CII 

As ubiquitous as Si II is in the spectra of SNe la, C II is 
nearly as rare. In a few cases, mainly at very early times, 
weak C II has been detected in optical spectra of SNe la 
(e.g.. Branch et al. 2003; Thomas et al. 2007; Tanaka et al. 
2008; Foley et al. 2010). Candidate SC SNe la, on the other 
hand, exhibit strong C II in their near-maximum spectra 
(Howell et al. 2006; Hicken et al. 2007; Scalzo et al. 2010). 

In our spectrum of SN 2009dc obtained 7 days before 
maximum brightness, we detect absorption from C II A6580 
and A7234 (both of which were also detected in a spectrum 
obtained the same day by Marion et al. 2009). The minimum 
of the absorption from C II A6580 is blueshifted by about 
9500 km s^^ in our pre-maximum spectrum. The expansion 
velocity of this feature is ~8500 km s~^ at 3 days before 

Note tliat they consider SNe la witli photosplieric velocities 
greater than 3iT above the average to be "high-velocity SNe la." 
^ This second velocity has a large uncertainty since it was mea- 
sured by approximating the minimum of the Si II absorption fea- 
ture by eye from Figure 3 of Contreras et al. (2010). 
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Figure 6. Our prc-maximum spectrum of SN 2009dc (and various comparison SNe) witli a few major lines identified and days relative 
to maximum light indicated for each spectrum (in parentheses). From top to bottom, the SNe are the "standard overluminous" Type la 
SN 1991T (Filippenko et al. 1992), a possible SC SN la SN 2004gu (Contreras et al. 2010, though this spectrum is from our own database), 
another possible SC SN la SN 2006gz (Hickcn et al. 2007) at two different epochs, SN 2009dc (this work), yet another possible SC SN la 
SN 2003fg (Howell et al. 2006), the SN 2002cx-like pecuHar SN 2005hk (Chornock et al. 2006), and finally the "standard normal" Type la 
SN 2005cf (Wang et al. 2009b). All spectra have had their host-galaxy recession velocities removed and have been dereddened according 
to the values presented in their respective references. Note that Si II, S II, and O I are often found in near-maximum spectra of SNe la, 
but C II is not. 



maximum, slowing to about 7000 km s ^ at 3 days after 
maximum, and ~6700 km by 6 days after maximum 
(Yamanaka et al. 2009; Tanaka et al. 2010). By about 18 
days past maximum this feature seems to have completely 
disappeared, though it is difficult to be sure given the low 
S/N spectrum seen in Figure 3 of Yamanaka et al. (2009). 
However, it is clear from the same figure that by 25 days 
after maximum, the feature is most definitely gone. It is also 
undetected in our spectrum obtained 35 days past maximum 
(see §3.4). 

Also, from our pre-maximum spectrum of SN 2009dc, 
we measure the minimum of the C II A7234 absorption to 
be ~9400 km s~^, which is in good agreement with the 
velocity of the C II A6580 feature that we calculate above. 
This feature is still apparent at about 6 days past maximum, 



near a velocity of ~6600 km s^^ (sec Fig. 1 of Tanaka et al. 
2010, though the feature is not marked), but it has also 
Ukely disappeared by 18 days past maximum (Yamanaka 
et al. 2009, Fig. 3) and is certainly gone by 35 days past 
maximum (see §3.4). 

For C 11, the velocities seen in SN 2009dc once again 
match those of SN 2003fg at ~2 days past maximum (How- 
ell et al. 2006, Fig. 3 and Supplementary Information), 
though only the C 11 A4267 line (at an expansion veloc- 
ity of ~8300 km s~^) is clearly detected. On the other hand, 
SN 2007if appears to have both the C II A6580 and A7234 
features in its spectrum from 5 days past maximum (Scalzo 
et al. 2010). 

A spectrum of SN 2006gz at 13 days before maximum 
shows both the C II A6580 and C II A7234 Unes with expan- 
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sion velocities near 15,500 km s^^. However, both features 
are effectively undetectable 4 days later (Hicken et al. 2007). 
We also note that neither line is seen in spectra of SN 2004gu 
at 5 days before maximum (sec Fig. 6, although this spec- 
trum has low S/N) or at 3 days after maximum (Contreras 
et al. 2010).'' 

For SN 2009dc, Marion et al. (2009) reported no de- 
tectable absorption from C II A4745 in their spectrum ob- 
tained 7 days before maximum. However, in our spectrum 

taken on the same day, wc do sec evidence for absorption 
from this transition, as well as absorption from C II A4267 
(see §3.3.4). Note that both of these features are slightly 
blended with absorption from other species, which is what 
most likely led Marion et al. (2009) to conclude that the 
features were not present. 

The EW of the C II A6580 and A7234 absorptions, as 
calculated from our spectrum taken 7 days before maximum, 
arc ~19 A and ~13 A, respectively. Hicken et al. (2007) re- 
port that the C II A6580 feature in SN 2006gz 14 days before 
maximum had an EW that was slightly higher than these 
values (25 A). Both of these SNe have much stronger C II 
lines than typical SNe la. For example, one of the strongest 
C II features ever seen in a normal SN la was in SN 2006D. 
where C II A6580 was detected with an EW of merely 7 A 
(Thomas et al. 2007). 

3.3.3 Ca II 

In our pre-maximum spectrum (as well as most of our post- 
maximum spectra; see §3.4) we detect narrow, as well as 
broad, absorption from Ca II H&K. The narrow components 
are at the redshift of SN 2009dc and its host galaxy. At all 
epochs these lines are just at our instrumental resolution 
limit; thus, calculating an accurate FWHM is nearly im- 
possible. In addition, the H&K lines occur in a part of the 
spectrum where numerous broad absorption features blend 
together, so defining a local continuum is difficult, exacer- 
bating the problem of measuring reliable line strengths. De- 
spite all this, we attempt to measure the strength of this 
narrow feature. We calculate that the Ca II H&K lines ap- 
pear to have FWHM ranging from about 100-300 km s~^. 
However, the line strengths are consistent with no change 
during these epochs, given our rather large uncertainties. 

Some SNe la exhibit narrow, time-variable Na I D ab- 
sorption lines with unchanging, narrow Ca II H&K ab- 
sorption (e.g., Simon et al. 2009). However, time-variable 
Ca II H&K absorption has not been seen. SN 2009dc does 
show weak, narrow absorption from Na I D, but it is far 
too weak to determine if variability in the line strengths ex- 
ists (see §3.5). We cannot claim to have detected narrow, 
time-variable Ca II H&K absorption in SN 2009dc, but such 
variability should be searched for in other, future SC SNe la. 

The narrow Ca II H&K may arise from calcium-rich in- 
terstellar material (ISM) in the host galaxy of SN 2009dc, 
UGC 10064, and in fact this feature is strong in a spec- 
trum of the core of the galaxy (Abazajian et al. 2009). 
However, since SN 2009dc is relatively far from the host's 
nucleus, this seems somewhat unlikely. The alternative is 

* The C II A7234 line, if present, would be at the red end of this 
spectrum. 



that the calcium-rich material which gives rise to the strong 
Ca II H&K line is in close proximity to the SN site. This 
means that SN 2009dc likely exploded in a region of calcium- 
rich ISM, or that the progenitor star of SN 2009dc itself had 
calcium-rich circumstellar material. However, which of these 
two possibilities best reflects the true situation and how this 
may affect the SN itself are beyond the scope of this paper. 

Another Ca II feature commonly seen in SNe la is the 
broad Ca II near-IR triplet. However, we see no hint of this 
line in our spectrum obtained 7 days before maximum. On 
the other hand, Tanaka et al. (2010) clearly detect this fea- 
ture near 8400 A in their spectropolarimetric observations 
taken ~6 days after maximum. Not only is the line appar- 
ent in their total-fiux spectrum, but it is the most polar- 
ized spectral feature they detect, with a total polarization 
of ~ 0.7%. This, along with the high polarization of Si II 
A6355 and the relatively low continuum polarization, led 
Tanaka et al. (2010) to conclude that SN 2009dc had a 
roughly spherical photosphere (as is the case for most nor- 
mal SNe la, e.g., Leonard et al. 2005; Chornock & Filippenko 
2008) with a somewhat clumpy distribution of intermediate- 
mass elements (IMEs). 

3.3.4 Other Species and the SYNOW Fit 

To determine which other species are present in our pre- 
maximum spectrum of SN 2009dc, we use the spectrum- 
synthesis code SYNOW (Fisher et al. 1997). SYNOW is 
a parameterized resonance-scattering code which allows for 
the adjustment of optical depths, temperatures, and veloci- 
ties in order to help identify spectral features seen in SNe. 

Before fitting, we deredden our pre-maximum spec- 
trum of SN 2009dc using E(B - V)mw = 0.070 mag, 
E{B — V)host = 0.1 mag (see §3.5 for more information 
on how we derive these values), Rv = 3.1, and the red- 
dening curve of Cardelli ct al. (1989). We also remove the 
host-galaxy recession velocity {cz = 6300 ± 300 km s"'^, as 
determined from the narrow Ca II H&K absorption in our 
spectrum) before comparing with the output of SYNOW. 

Our derived SYNOW fit is compared to our actual pre- 
maximum spectrum in Figure 7. Also shown arc the spectral 
features from each of the individual species that were used in 
the final fit. Our SYNOW fit faithfully reproduces the vast 
majority of the features seen in our pre-maximum spectrum 
of SN 2009dc, but a few features are not matched exactly 
(specifically, the ones near 3700 A, 4100 A, and 4600 A). 

In addition to the species mentioned above (Si II, C II, 
and weak Ca II), we detect O I, S II, Mg II, Si III, Si IV, 
Fe II, Co II, and Co III. All of these ions (except C II, as 
mentioned previously) are often found in early-time spectra 
of SNe la, though O I may be weak in some overluminous 
SNe la (e.g., Fihppenko 1997). The detection of C II in the 
spectra of SN 2009dc (see §3.3.2) implies that it had a signif- 
icant amount of unburncd material in its cjccta, and thus the 
oxygen seen in SN 2009dc is most likely from this unburned 
material as well. 

Our SYNOW fit has a photospheric velocity of 
9000 km s"\ with maximum velocities of 10,000— 
15,000 km s"'^ depending on the ion. This value matches 
our derived velocities from the Si II and C II absorption fea- 
tures (see §3.3.1 and §3.3.2, respectively). Since the early- 
time spectra of SN 2009dc and SN 2003fg greatly resemble 
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Figure 7. The top cloven speetra show the eonstitucnt species in our SYNOW fit. The second from bottom spectrum is the sum of 
the top eleven spectra. The bottom spectrum is our pre-maximum spectrum of SN 2009dc (dereddcncd by E{B — y)MW = 0.070 mag 
and E(B — V)host =0.1 mag with Ry =3.1). We also remove the host-galaxy recession velocity (cz = 6300 ± 300 km s~^). Finally, we 
remove the underlying 20,000 K blackbody continuum from all spectra shown. 



each other, it is not surprising that these parameters are sim- 
ilar to those found in a SYNOW fit of SN 2003fg (Howell 
et al. 2006, Supplementary Information). Our fit requires a 
slightly higher photospheric velocity (they use 8000 km s"'^), 
but this seems reasonable since our spectrum was obtained 
7 days before maximum and the spectrum of SN 2003fg fit 
by Howell et al. (2006) is from ~2 days after max;imum. 

All ions in our SYNOW fit had an excitation temper- 
ature around 10,000 K except C II, which was raised to 
20,000 K in order to get the relative line strengths to match 
the data. This is lower than what was calculated by How- 
ell et al. (2006) for SN 2003fg. They required an excitation 
temperature for C II of 35,000 K to match their spectrum, 
while we obtain a temperature that is less than 60% of that 
value. 

Finally, we require an underlying blackbody tempera- 
ture of 20,000 K in order to reproduce the overall continuum 
shape of SN 2009dc. This is hot for a SN la, even at early 



times, and is indicative of the production of a large amount 
of ^''Ni (Nugent et al. 1995). Blackbody temperatures in 
the range of 10,000-15,000 K seem to be more common for 
normal SNe la (e.g., Patat et al. 1996), with the overlumi- 
nous Type la SN 1991T having temperatures at the high 
end of that range (Mazzafi et al. 1995). Howell et al. (2006) 
used a blackbody temperature of 9000 K to fit SN 2003fg, 
which does not seem abnormally high for a SN la 2 days 
past majcimum. Even though this value is a bit lower than 
our derived blackbody temperature, it is not too surprising 
since our spectrum of SN 2009dc was obtained about 9 days 
earlier (relative to maximum light) than their spectrum of 
SN 2003fg. 



3.4 Post-Maximum Spectra 

The rest of our spectral data of SN 2009dc (i.e., all of our 
post-maximum spectra) are presented in Figure 8. Note that 
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during our observation on 2009 Aug. 14.2 (day 109), there 
was smoke and ash in the sky throughout the night (from a 
fire near the observatory), and thus the continuum shape is 
less accurate (due to the nonstandard extinction caused by 
the smoke) than in the other observations. 

3.4.1 1 Month After Maximum 

Our day 35 (2009 May 31.3) spectrum of SN 2009dc is shown 
in Figure 9, along with spectra of the peculiar SN 2002cx 
(Li et al. 2003), another possible SC SN la SN 2007if (Scalzo 
et al. 2010, though the displayed spectrum is from our own 
database), and the "standard normal" Type la SN 2005cf 
(Wang et al. 2009b). While there are some spectroscopic 
similarities between SN 2009dc and SN 2005cf at this epoch, 
there are obvious differences as well. Most notably, there 
are many more medium-width features in the spectrum of 
SN 2009dc compared with SN 2005cf, and the features that 
do clearly match are at much lower velocities in SN 2009dc. 
This is the same as what was seen in the pre-maximum 
spectra (§3.3). Furthermore, the Ca II H&K feature as well 
as the Ca II near-IR triplet are much stronger in SN 2005cf 
than they are in SN 2009dc, again consistent with the pre- 
maximum spectrum of SN 2009dc (§3.3.3). 

The possible SC SN la SN 2007if appears to be a rea- 
sonable match to the spectrum of SN 2009dc, but there are 
some significant differences. Even though weak Ca II fea- 
tures are detected in both of these objects, there arc still a 
few medium-width features in SN 2009dc that are not seen 
in SN 2007if. In addition, the velocities of some (but not 
all) of the features in SN 2007if arc larger than those in 
SN 2009dc by ~1000 km s"\ However, SN 2007ff appears 
to match SN 2005cf well at this epoch, even though its ex- 
pansion velocity is still not as large as that of SN 2005cf, 
and SN 2007if is missing the Si II A6355 absorption seen in 
SN 2005cf. 

The best match to SN 2009dc at this epoch appears to 
be the peculiar SN 2002cx at an age about 5-10 days younger 
than SN 2009dc. A detailed discussion of the comparison be- 
tween peculiar SN 2002cx-like SNe and SN 2009dc can be 
found in §3.4.5. Since nearly every single spectral feature of 
SN 2009dc at this epoch matches those in SN 2002cx, we use 
the line identifications of Li et al. (2003) (which mainly came 
originally from Mazzah et al. 1997) for SN 2009dc. The spec- 
tra at this time are dominated by various multiplcts of Fc II, 
although Ca II (as mentioned above), Co II, and Na I are 
also present. There is evidence for O I A7773 and O I A9264 
(which may give rise to the features seen near 7600 A and 
9100 A, respectively), though Branch et al. (2004) claim that 
a mixture of Fc II and Co II transitions can account for these 
parts of the spectrum without the need for oxygen. Possible 
detections of Ti II (Li et al. 2003) and Cr II (Branch et al. 
2004) have been claimed for SN 2002cx, but it is extremely 
difficult to disentangle the individual contributions of each 
of these IGEs in SN 2009dc. 

The one clear difference between the spectrum of 
SN 2009dc and that of SN 2002cx (and SN 2005hk as well) 
is that SN 2009dc has a strong absorption feature centred 
at 6250 A that is not seen in cither of these other pecu- 
liar SNe la. This line is also not present in our spectrum 
of SN 2007if from a similar epoch. In Figure 10 we present 
a SYNOW fit to our spectrum of SN 2009dc from 35 days 



past maximum based originally on a fit to the spectrum of 
SN 2002cx from 25 days past maximum (Li et al. 2010). 
Both SYNOW fits reproduce the majority of the strongest 
features in SN 2009dc, and the only difference between the 
two is the addition of Si II at an expansion velocity and ex- 
citation temperature similar to other ions in the fit. Based 
on the comparison of our two fits, we attribute this myste- 
rious absorption in SN 2009dc to Si II A6355 (blueshifted 
by ~5500 km s"^). Note that the spectrum of SN 2005cf 
from 28 days past maximum, shown in Figure 9, has strong 
Si II A6355 absorption as well. The main difference between 
these two features is that the absorption in SN 2005cf is 
blueshifted by ~9100 km s~^, nearly twice the velocity seen 
in SN 2009dc. We discuss this feature further in §3.4.3. 



3.4.2 2 Months After Maximum 

Our day 52 (2009 Jun. 17.5) and day 64 (2009 Jun. 29.3) 

spectra of SN 2009dc are shown in Figure 11, along with 
spectra of SN 2002cx (Li et al. 2003) and SN 2007ff (Scalzo 
et al. 2010, but the displayed spectrum is from our own 
database) . Note that SN 2009dc does not evolve much spec- 
troscopically between days 52 and 64, though there are a 
few differences which will be discussed in §3.4.3. 

As was the case a month earlier, the spectrum of 
SN 2009dc at 2 months past maximum matches that of 
SN 2002cx at a similar epoch. Again, the overall shape and 
line profiles look extremely similar between these two objects 
for the most part, and both look similar to SN 2007if as well 
as SN 2005hk (Phillips et al. 2007). At this epoch there re- 
main a handful of medium- width features in SN 2009dc that 
are not seen in SN 2007if, and once again the velocities of 
some features in SN 2007if appear to be larger than those 
of SN 2009dc by about 1000 km s~\ 

As mentioned above, the Si II A6355 absorption that 
was seen in our day 35 spectrum is still present in our day 
52 and day 64 spectra of SN 2009dc. However, it decreases 
in strength significantly during these three observations and 
disappears almost completely by 80 days past maxinnim (see 
§3.4.3). Another major difference between SN 2009dc and 
SN 2002cx is that even though the Ca II H&K features are 
similar in strength, emission from the Ca II near-IR triplet 
is weaker in SN 2009dc and has a complex profile with mul- 
tiple peaks. We discuss this feature further in the following 
section. 



3.4.3 3-4 Months After Maximum 

Figure 8 shows all of our post-maximum spectra of 
SN 2009dc. Overall, the spectra do not change very much 
from 35 to 109 days past maximum (see §3.4.4 for more infor- 
mation on our spectrum obtained 281 days past maximum). 
Most of the major spectral features remain constant during 
these epochs (except for a few notable examples discussed 
below) and their profile shapes are largely unaltered as well. 
The profiles of the features that do change with time mainly 
seem to weaken in a symmetric manner. That is, there is no 
compelling evidence for a preferential increase or decrease 
of fiux in the blue or red wings of any feature. 

One of the most striking aspects of the spectral evolu- 
tion of SN 2009dc over these epochs is the disappearance 
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Figure 8. Our post-maximum spectra of SN 2009dc, witii days relative to maximum liglit indicated for each spectrum. Important 
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of the broad emission near 5400 A and in the ranges 6500- 
7200 A, 7800-8300 A, and (to a lesser extent) 8800-9200 A. 
The first two of the ranges mainly fall into the R and / 
bands, respectively. The emission from these features, which 
is strong through about 1 month past maximum and then 
begins to decrease with time, is possibly responsible for the 
platcau-likc broad peak in the R and I bands near maximum 
brightness and their slow late-time decline (as compared to 
more normal SNe la, see Fig. 2). 

Li et al. (2003) also saw these features in SN 2002cx, 
and their temporal evolution and the effect they had on the 
near-maximum R and 7-band light curves is nearly identical 
to what we find in SN 2009dc. However, Li et al. (2003) were 
unable to determine which ions were responsible for some of 
the features in these ranges in SN 2002cx, and we are like- 
wise unsure of their origin. The detailed spectral modeling 
required to definitively determine their identity is beyond 
the scope of this paper. 

As noted previously, we detect strong Si II A6355 ab- 



sorption in our day 35 spectrum of SN 2009dc. This feature 
is still present, but weakens, in our spectra from 52 and 
64 days past maximum. In our observation of SN 2009dc at 
80 days past maximum there remains a hint of absorption 
from Si II A6355, but it is almost nonexistent and it is effec- 
tively undetected in any of our later spectra. It is unusual 
to sec such strong absorption from Si II A6355 in a SN la at 
such late times. For example, at 46 days past majcimum the 
spectrum of SN 2005cf has only a barely perceptible Si II 
A6355 feature (Wang et al. 2009b, Fig. 17). 

The relative strength of this Si II line at late times in 
SN 2009dc is possibly due to its large cjecta mass. In §4.2 
we find that SN 2009dc has one of the largest ^^Ni masses 
ever produced by a SN la, and thus one of the largest total 
ejecta masses as well. The massive ejecta may imply that 
a greater than average amount of silicon was synthesised 
in the explosion of SN 2009dc, which could give rise to the 
rather long-lived Si II A6355 line that we detect. However, it 
is also possible that the Si-rich ejecta have an aspherical or 
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Figure 9. Our spectrum of SN 2009dc 35 days past maximum and a few comparison SNe, with days relative to maximum liglit indicated 
for each spectrum (in parentheses). From top to bottom, the SNe are the peculiar SN 2002cx (Li et al. 2003, where we have combined 
their spectra from 20 and 25 days past maximum), SN 2009dc (this work), another possible SC SN la SN 2007if (Scalzo et al. 2010, 
though the spectrmn shown is from our own database), and the "standard normal" Type la SN 2005cf (Wang et al. 2009b). Important 
features discussed in the text are labeled. All spectra have had their host-galaxy recession velocities removed and have been dereddened 
according to the values presented in their respective references. 



clumpy distribution, and thus the longevity of the observed 
Si II A6355 line is simply a viewing angle effect. 

As stated earlier, another major difference between 
SN 2009dc and SN 2002cx is that emission from the Ca II 
near-IR triplet is weaker in SN 2009dc and has a complex 
profile with multiple peaks, even though the Ca II H&K fea- 
tures are similar in strength. A SYNOW fit indicates that we 
are actually resolving the individual components of the Ca II 
near-IR triplet along with shallow absorption possibly from 
Mg I and another, unidentified ion. The Ca II and Mg I fea- 
tures are blueshifted by ~6000 km s~^, which matches fairly 
well with our previously calculated post-maximum expan- 
sion velocities (see §3.4.2). Thus, it seems that the complex 
structure seen in our post-maximum spectra of SN 2009dc 
near 8600 A is a mixture of the Ca II near-IR triplet, Mg I, 
and probably another ion whose identity we are unable to 
definitively determine. 



3.4-4 9 Months After Maximum 

On 2010 Feb. 6.6 (281 days after maximum), we obtained 
our final spectrum of SN 2009dc, shown in Figure 12 along 
with late-time spectra of the normal Type la SN 1998bu (Jha 
et al. 2006), SN 2006gz (Maeda et al. 2009), another normal 
Type la SN 1990N (Gomez & Lopez 1998, via the online 
SUSPECT database"'), and SN 2002cx (Jha et al. 2006). 

At over 9 months past maximum light, SN 2009dc ap- 
pears to have completely transitioned to the nebular phase. 
The SN ejecta are now optically thin and the spectrum 
shows various emission features. The broad peaks in the 
range 3800-5500 A arc usually attributed to blends of var- 
ious features from [Fe II] and [Fe III] (e.g., Mazzali et al. 
1998). These blends appear slightly weaker in SN 2009dc as 
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Figure 10. The top panel shows our spectrum of SN 2009dc from 35 days past maximum (solid line) overplotted witli our best SYNOW 
fit which includes Si II (dotted line). The bottom panel is the same as the top panel but zoomed in and centred on the absorption feature 
in SN 2009dc near 6250 A; it also shows the same SYNOW fit but with no Si II (dashed line). The SYNOW fit with Si II appears to 
reproduce the data better, and thus we attribute this absorption to Si II A6355 (blueshifted by ^5500 km s~^). We have dereddened 
our data by E(B — V)mw = 0.070 mag and E(B — y)host = 0-1 mag with Ry = 3.1. We have also removed the host-galaxy recession 
velocity (cz = 6300 ± 300 km s~^) from the data. 



compared to the normal Type la SNc 1998bu and 1990N 
at similar epochs. Interestingly, Maeda et al. (2009) found 
that SN 2006gz completely lacks these features, though their 
spectrum has a low S/N. 

The Ca II near-IR triplet is detected in both of the nor- 
mal SNe la presented in Figure 12 as well as in SN 2009(lc 
and possibly SN 2006gz (though the spectrum is noisy at 
these wavelengths). There is essentially no emission in ei- 
ther SN 2009dc or SN 2006gz from ~5700 A to ~6800 A, 
a region in which [Co III] emission has been identified in 
normal SNe la (Kuchner et al. 1994; Ruiz-Lapuente et al. 
1995). The ratio of emission from the strongest [Fe III] fea- 
ture (4589-4805 A) to emission from the strongest [Co III] 
feature (5801 5995 A) has been measured in a handful of 
SNe la at late times; it increases with both phase relative to 
maximum and ejecta temperature (Kuchner et al. 1994). 

Using our spectra of SNe 1998bu and 1990N shown in 



Figure 12, we calculate ratios of about 5.6 and 7.7, respec- 
tively. These fit nicely with the other data displayed in Fig- 
ure 2 of Kuchner et al. (1994). However, it is nearly impos- 
sible to calculate this ratio for SN 2009dc and SN 2006gz, 
since they appear to lack almost any putative cobalt emis- 
sion at this epoch (and the cobalt and iron lines in spectra 
at earlier epochs are still significantly blended). We attempt 
to measure the tiny amount of [Co III] emission above the 
background level from 5801 A to 5995 A and calculate a 
Fe/Co ratio of about 27. This value, which is likely an un- 
derestimate, is already well above that of any other SN la 
at this phase presented by Kuchner et al. (1994). It is also 
above the dotted lines in Figure 2 of Kuchner et al. (1994), 
implying that the ejecta temperature of SN 2009dc may be 
> 10, 000 K. Even though this is broadly consistent with our 
relatively largo value of 20,000 K for the blackbody temper- 
ature of SN 2009dc before maximum (§3.3.4), 10,000 K is 
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Figure 11. Our spectra of SN 2009dc about 52 and 64 days past maxiiiiuiii, and a few comparison SNe with days relative to maximum 
light indicated for each spectrum (in parentheses). The top SN is tiie peculiar SN 2002cx (Li ct al. 2003), the next two spectra are 
SN 2009dc (this work), and the bottom spectrum is another possible SC SN la SN 2007if (Scalzo et al. 2010, though the displayed 
spectrum is from our own dataliase). Important features discussed in the text are labeled. All spectra have had their host-galajcy 
recession velocities removed and have been dereddened according to the values presented in their respective references. 



extremely hot for a normal SN la 9 months past maximum 
brightness and the shape of the spectrum at this epoch does 
not seem to support such a high temperature. 

The most well-defined spectral feature in our day 281 
spectrum of SN 2009dc is the broad, doublc-pcaked line cen- 
tred near 7200 A, which was also the strongest feature de- 
tected in a spectrum of SN 2006gz from 338 days past max- 
imum (Maeda et al. 2009). In late-time spectra of SNe la 
this is attributed to blends of various [Fe II] lines, whereas in 
coro-coUapso SNe it is associated with [Ca II] AA7291, 7324, 
though in a lato-timc spectrum of SN 2006gz it has been 
suggested that the feature is due to a combination of both 
[Fe II] and [Ca II] hues (Maeda et al. 2009). The shape of 
the profile in SN 2009dc matches roughly that of SN 1998bu 
(though the central dip in SN 2009dc is deeper), but differs 
somewhat significantly from that of SN 1990N. Shapes simi- 
lar to the one seen in SN 2009dc have, however, been seen by 
Maeda et al. (2009), who present synthetic late-time spectra 



of SNe la. Their models with more massive progenitors seem 
to lead to the feature near 7200 A becoming strong relative 
to the complex of lines in the range 3800-5500 A, as seen in 
SN 2009dc. 



3.4.5 SN 2009dc and SN 2002cx-Like Objects 

As already mentioned numerous times, SN 2009dc shares 
some interesting properties with both SN 2002cx and 
SN 2005hk, which arc members of a class of peculiar SNe la 
(Fihppenko 2003; Li et al. 2003; Jha et al. 2006). These so- 
called "SN 2002cx-like" SNe are characterized by low expan- 
sion velocities and low peak luminosities, slow late-time de- 
clines in R, and Fe III features dominating early-time spectra 
(Jha et al. 2006). Many of these criteria hold for SN 2009dc 
as well, but the very obvious difference is that the peak lu- 
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Figure 12. Our spectrum of SN 2009dc about 281 days past maximum, and a few eomparison SNe with days relative to maximum light 
indicated for each spectrum (in parentheses). The spectra from top to bottom arc the normal Type la SN 1998bu (Jha et al. 2006), 
SN 2009dc, SN 2006gz (Maeda et al. 2009), SN 1990N (via SUSPECT), and the pecuUar Type la SN 2002cx (Jha et al. 2006). Major 
spectral features are labeled. All spectra have had their host-galaxy recession velocities removed and have been dereddened according to 
the values presented in their respective references. 



minosity of SN 2009dc is much larger than average (see §4.1) 
while those of SN 2002cx-likc objects arc well below average. 

In §3.1, we noted the absence of a prominent second 
maximum in the R and /-band light curves of SN 2009dc. 
This lack of a secondary maximum has boon observed in 
SN 2002cx-like objects (Li et al. 2003; Phillips ct al. 2007). 
However, as seen in Figure 2, SN 2005hk (a SN 2002cx-like 
object) evolves much faster than SN 2009dc in all bands 
and proves to be a poor match photometrically. The colour 
indices of SN 2009dc and SN 2005hk are again a poor match 
as shown in Figure 5. 

Spcctroscopically, SN 2009dc does not really resemble 
SN 2002cx-like objects before maximum (see Fig. 6), nor 
about 9 months after maximum (sec Fig. 12). At these lat- 
est times, the spectrum of SN 2009dc consists of broad 
emission from forbidden iron lines while the spectrum of 
SN 2002cx is made up of many narrow permitted lines of 
Fe II, Na II, Ca II, and (tentatively) O I (Jha et al. 2006). 



This is supported by the fact that the narrow spectral fea- 
tures of SN 2002cx do not simply appear to correspond to 
resolved versions of the blends seen in SN 2009dc. 

Despite all of these differences, we have shown that 
at 1-2 months past maximum brightness the spectra of 
SN 2009dc and SN 2002cx arc surprisingly similar. In Fig- 
ure 9 we have combined the day 20 and 25 spectra of 
SN 2002cx from Li et al. (2003) in order to improve the 
S/N and the wavelength coverage. Not only are the overall 
shapes of these two spectra nearly identical, but the width, 
strength, and existence of almost all of the individual spec- 
tral features match extremely well (with one significant ex- 
ception, the Si II A6355 absorption as mentioned in §3.4.1). 
Both objects, after being dereddened and deredshifted, ap- 
pear to have an underlying blackbody continuum of about 
5500 K (Branch et al. 2004). 

About a month past maximum, Li et al. (2003) and 
Branch et al. (2004) calculate expansion velocities of ~5000- 



© 2010 RAS, MNRAS 000, 1-30 



SN 2009dc, Super-Chandrasekhar Mass? 19 



6000 km s"^ for SN 2002cx, and Phillips et al. (2007) de- 
termine expansion velocities in nearly the same range for 
SN 2005hk. For SN 2009dc we also derive velocities in this 
range for a few ions (O I, Ca II, Fe II, and Si II) with sig- 
nificant uncertainty since not only arc many of those fea- 
tures blended, but the line identifications themselves are not 
definitive. These values also match well with the expansion 
velocity of ~6000 km derived by Yamanaka et al. (2009) 
from their spectrum of SN 2009dc from 25 days past maxi- 
mum. 

Both Li et al. (2003) and Branch et al. (2004) deter- 
mine an expansion velocity for SN 2002cx at 56 days past 
maximum of ~4700 km s~^ based on the Fc II A4555 feature 
and a SYNOW fit. Phillips et al. (2007) calculate a slightly 
higher value of about 6000 km s"^ for SN 2005hk at 55 days 
past maximum (also using the Fe II A4555 line). Based on 
this Fe II feature, we determine an expansion velocity of 
~5500 km s-i for our spectra of SN 2009dc from 52 and 
64 days past maximum. If we use the Si II A6355 line, which 
is still present in SN 2009dc but not seen in the SN 2002cx- 
like objects, we derive similar expansion velocities for both 
our day 52 and day 64 spectra. In general, it appears that 
the post-maximum expansion velocities of the spectral fea- 
tures of SN 2009dc are slightly greater than those seen in 
SN 2002cx by a few hundred km s~^ depending on the fea- 
ture (Li et al. 2003; Branch et al. 2004), whereas they are 
about equal to those of SN 2005hk (Philhps et al. 2007). 

In conclusion, SN 2009dc and SN 2002cx-likc objects 
have quite different photometric properties but are nearly 
identical spectroscopically only at 1-2 months past maxi- 
mum brightness. The expansion velocity, the blackbody tem- 
perature, and the chemical composition match extremely 
well at those epochs, but the evolution of these observablcs 
from pre-maximum to nearly a year after maximum is quite 
diflferent between SN 2009dc and SN 2002cx-Uke objects. 

3.4.6 SN 2007if 

In Figure 13 we present our late-time spectra of another 
possible SC SN la, SN 2007if. As noted earlier, the spectra 
of SN 2007if generally look similar to those of SN 2009dc, 
SN 2002cx, and even SN 2005cf at comparable epochs. The 
most significant spectral differences between SN 2007if and 
SN 2009dc arc, as stated previously, that the expansion ve- 
locities seen in SN 2007if arc '^lOOO km s^^ larger than those 
in SN 2009dc and SN 2002cx at comparable phases. This 
larger velocity is likely causing some of the medium-width 
features seen in SN 2009dc and SN 2002cx to be "smoothed 
out" , and hence we do not detect many of these features in 
SN 2007if. It also has velocities which arc smaller than those 
seen in SN 2005cf, but otherwise the spectra match well. 

We have also stated previously that SN 2007if does not 
show the strong Si II A6355 absorption after maximum that 
is visible in spectra of SN 2009dc (as well as in SN 2005cf). 
Unfortunately, our observations of SN 2007if do not extend 
to sufficiently long wavelengths for us to inspect whether 
the Ca II near-IR triplet of SN 2007if shows the complex 
structure seen in SN 2009dc. 

The main aspects of the spectral evolution of SN 2007if 
over our three epochs is the disappearance of broad emis- 
sion near 5500 A and in the ranges 6700-7100 A and possi- 
bly 7800-8100 A (though our data are noisy in this second 



range of wavelengths). These match nearly exactly with the 
spectral features of SN 2009dc that were seen to decrease in 
strength over similar epochs.^ 

3.5 Reddening and Extinction 

As mentioned above, we detect weak (but noticeable), nar- 
row absorption from Na I D in all of our spectra of 
SN 2009dc, both at zero redshift and at the recession ve- 
locity of the SN itself {cz = 6300 ± 300 km s"^). Even 
though there arc clear notches in all of our spectra at the 
correct wavelengths for Na I D in the MW and in the host 
of SN 2009dc (UGC 10064), they are not broader than our 
spectral resolution. 

Wo attempt to measure an EW of both absorption fea- 
tures in each of our observations, and we calculate average 
values of ~0.7 ± 0.1 A for the MW and ~1.0 ± 0.1 A for 
UGC 10064. These are close to the values of 0.7 A and 1.2 A 
reported by Marion et al. (2009) and 0.5 A and 1.0 A found 
by Tanaka et al. (2010). 

There is some uncertainty in how exactly to convert the 
EW of Na I D absorption into extinction, but a few (mainly 
empirical) relations do exist (e.g., Barbon et al. 1990; Mu- 
nari & Zwitter 1997; Turatto et al. 2003). Using these con- 
versions, we determine a range of values for E{B — V) both 
in the MW (0.1 < E{B -V)mw < 0.18 mag) and in 
UGC 10064 (0.15 < E{B - l/)ho.t < 0.3 mag). 

For the extinction from the MW, we compare our cal- 
culated values of E(B — V) to the value from the Galac- 
tic dust maps of Schlegel et al. (1998). Their maps give 
E{B - 1/)mw = 0.070 mag for the position of SN 2009dc, 
slightly lower than our range of values. However, since the 
accuracy of the Schlegel et al. (1998) dust maps is almost cer- 
tainly higher than that of our EW measurements, we adopt 
E{B — 1^)mw = 0.070 mag for our analysis. 

The extinction from UGC 10064 is a bit more com- 
plicated to determine accurately. If the Lira-Phillips rela- 
tion holds for SN 2009dc, then we expect an extinction of 
E{B — V\ost ~ 0.26 mag (see §3.2 for details on our ap- 
plication of the Lira-Phillips relation), which is within our 
range of values calculated from the spectra. Furthermore, 
if we assume that the EW of Na I D is directly propor- 
tional to E{B — V), as Tanaka et al. (2010) do, then we 
derive E{B - V^)host ~ (l.O A/0.7 A) x 0.070 « 0.1 mag. 
Taking into account all of these values, their uncertain- 
ties, and the method by which they were derived, we adopt 
E{B — V^)host = 0.1 mag for our analysis, which is our low- 
est plausible, nonzero reddening. However, we note that the 
reddening could be as large as 0.3 mag. 

3.6 Host-Galaxy Properties 

SN 2009dc occurred ~27" away from the centre of the SO 
galaxy UGC 10064. We adopt the luminosity distance {dh) 
as the distance to the SN, again assuming a ACDM Universe 
with Ho = 70 km s"^ Mpc"\ Qm = 0.27, and Qa = 0.73 



^ The third wavelength range which possibly decreased in 
strength in SN 2009dc, 7800-8100 A, is outside of our spectral 
coverage of SN 2007if. 
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Figure 13. Our spectra of SN 2007if with days relative to maximum lighit indicated in each case. We have removed the host-, 
recession velocity (cz = 22,200 ±250 km s"!, Scalzo et al. 2010). 



^iaxy 



(Spergel et al. 2007) and z = 0.022 ± 0.001, where the red- 
shift has been corrected for infall into the Virgo cluster and 
we have adopted an uncertainty of cz = 300 km s^^ to 
allow for a peculiar velocity induced by the gravitational 
interaction of neighboring galaxies. At di = 97 ± 7 Mpc, 
SN 2009dc is a projected distance of ~ 13 kpc from the 
centre of UGC 10064. A useful galactic scale is the ra- 
dius containing 50% of the flux (Petrosian 1976), which for 
UGC 10064 is about 4'.'71 ± 0'.'21 in the R band (Abaza- 
jian et al. 2009). This puts SN 2009dc about 5.7 Petrosian 
radii (in projection) from the centre of its host. For compar- 
ison, SN 2003fg was found to be projected 0.9 kpc from the 
core of its low-mass star-forming host (Howell et al. 2006) 
and SN 2007if was discovered effectively directly on top of 
its low-mass and low-metallicity host (Scalzo et al. 2010). 
The host of SN 2006gz is an Scd galaxy, and the SN was 
projected 14.4 kpc from the nucleus (Hicken et al. 2007), 
while SN 2004gu was projected 2.21 kpc from its host's nu- 
cleus (Monard et al. 2004). It seems that host-galaxy type 
and distance from the nucleus vary widely among candidate 



SC SN la (though the distance comparison is difficult given 
projection effects). However, SN 2003fg, SN 2007if, and pos- 
sibly SN 2009dc (see below) were all associated with regions 
of recent star formation. 

Wegner & Grogin (2008) inspected an optical spectrum 
of the core of UGC 10064 and derived an age of 4.0±3.0 Gyr, 
a metallicity [Z /H] = 0.45±0.07, and a chemical abundance 
[a/Fe] = 0.13 ± 0.09. This value for metallicity is about la 
above the average found for 26 SO/E void galaxies by Wegner 
& Grogin (2008), whereas the chemical abundance is just 
about the same as their average value. Optical spectra of 
the nucleus of UGC 10064 show no strong emission lines 
(Wegner & Grogin 2008; Abazajian et al. 2009), implying 
very little ongoing star formation (< 0.2 M© yr~^) in the 
core of the galaxy. 

However, recent deep H I surveys of SO galaxies have 
shown that the majority of nonclustcr early-type galaxies 
contain substantial reservoirs (lO'^-lO^ Mq) of neutral gas 
that extend to large galactocentric radii (Morganti et al. 
2006), whose morphology and kinematics suggest external 
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origins. On the other hand, in the majority of these galaxies 
the column density of this gas is too low to support star 
formation at large radii, and thus the stellar populations 
throughout the galaxy are old. 

Interestingly, in a fraction of the early-type galaxies rich 
in H I, a dynamical feature such as a prominent bar can cre- 
ate relatively dense ring-like accumulations of neutral gas at 
large radii (e.g., Schirmcrcr & Scovillo 2002; Weijmans et al. 
2008). In these regions of enhanced H I, low-level star forma- 
tion (~0.1 Mq yr~^) is observed (Jeong et al. 2007; Shapiro 
et al. 2010). Since SN 2009dc exploded in the outskirts of 
an SO galaxy, it seems difficult to say with certainty whether 
SN 2009dc came from an old or young stellar population. 

In the Supplementary Information of Howell et al. 

(2006), the age of the host galaxy of SN 2003fg is esti- 
mated to be ~700 Myr (with significant uncertainty) and 
the star-formation rate to be 1.26t.1'^^ yr '^i averaged 
over 0.5 Gyr. This star-formation rate is nearly an order of 
magnitude greater that the value derived by Wegner & Gro- 
gin (2008) for the core of UGC 10064 as well as the residual 
star-formation rate seen in rings of enhanced neutral gas at 
large galactocentric radii (Jeong et al. 2007; Shapiro et al. 
2010). SC SNe la are expected to be found preferentially in 
young stellar populations (e.g., Howell et al. 2006; Chen & 
Li 2009, and references therein), but they are not necessarily 
restricted to them. If SC SNe la are the result of the merger 
of two degenerate objects, then the merger timescale can be 
set by the emission of gravitational waves (e.g., Iben & Tu- 
tukov 1984) . Depending on the particular orbital parameters 
of the binary, there may be a significant delay between the 
formation of the constituent degenerate objects and their 
eventual merger. Thus, SC SNe la could occur in older stel- 
lar populations. 

Even though there is a significant difference between 
the average star-formation rates near the sites of SN 2009dc 
and SN 2003fg, there is still a possibility that they both 
come from relatively young stellar populations. As men- 
tioned above, an age of ~700 Myr was calculated for the 
environment of SN 2003fg (Howell et al. 2006, Supplemen- 
tary Information), while an age of ~4.0 Gyr was calculated 
for the core of the host of SN 2009dc (Wegner & Grogin 
2008). The fact that SN 2009dc is a few effective radii away 
from the nucleus of its host, coupled with the observations 
of low levels of residual star formation in other early-type 
galaxies at such galactocentric radii, imply that it is possi- 
ble that SN 2009dc also came from a relatively small, but 
young, stellar population on the edge of UGC 10064, and 
that the age calculated for the core of the galaxy is not ac- 
tually representative of the local environment of SN 2009dc. 

Furthermore, in the SDSS Data Release 7, an ex- 
tended, blue source appears 99'.'3 (which is a projected 
distance of 43.6 kpc) to the northwest of UGC 10064 
(Abazajian ct al. 2009). This SBdm galaxy is known as 
SDSS J155108.42+254321.3 or UGC 10063 (de Vaucoulcurs 
et al. 1991). Swift/XJVOT images of the SN 2009dc field 
(shown in Fig. 14) indicate that UGC 10063 is significantly 
bluer than UGC 10064. The irregularly shaped UGC 10063 
is barely visible in the KAIT _R-band image to the northwest 
of SN 2009dc in Figure 1. Stacked UVOT images show that 
UGC 10063 increases in prominence in progressively bluer 
bands. As seen in Figure 14, UGC 10063 seems to have a 
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Figure 14. Keck g-band image and stacked UVOT U, uvm2, 
and uvw2 images of UGC 10063 (SDSS J155108.42-I-254321.3). 
UGC 10063 and SN 2009dc are labeled in the upper-left and 
lower-right panels. The uijpcr-lcft panel is 2.69' X 2.69', and the 
remaining panels arc 4.67'x4.67 '; north is up and cast is to the 
left in all four images. UGC 10063 increases in prominence in 
progressively bluer bands, and is much bluer than UGC 10064 
(compare its appearance here with that in Fig. 1). It also seems 
to have a tidal tail extending toward UGC 10064, pointed roughly 
in the direction of the site of SN 2009dc. 



tidal tail extending toward UGC 10064, pointed roughly at 
the site of SN 2009dc. 

On 2010 Feb. 7 we obtained a spectrum of UGC 10063 
using LRIS on the Keck I telescope (with a l"-wide sUt, a 
600/4000 grism on the blue side, and a 400/8500 grating 
on the red side, resulting in FWHM resolutions of ~4 A 
and ~6 A, respectively). It has a blue continuum with few 
spectral features, but there is strong and narrow absorption 
from the Balmer series (though Ha is undetected). Strong 
Balmer lines and the lack of obvious emission lines are the 
hallmarks of so-called "post-starburst" galaxies in which vig- 
orous star- formation activity ended between ~50 Myr and 
1.5 Gyr ago (Poggianti et al. 2009, and references therein). 
This is consistent with UGC 10063 being a SBdm galaxy 
(de Vaucouleurs et al. 1991). From the seven Balmer ab- 
sorption lines marked in Figure 15, we calculate the redshift 
of UGC 10063 to be 0.021 ± 0.001, consistent with the red- 
shift derived from the 21-cm H I line {z = 0.02158±0.00003, 
Schneider et al. 1990). This also exactly matches the redshift 
of UGC 10064 and SN 2009dc (as determined from narrow 
Ca II H&K absorption). 

Based on the spectrum of UGC 10063 and its position 
relative to UGC 10064, we propose that these two galax- 
ies had a close encounter in the relatively recent past. If 
UGC 10063 passed near the nucleus of UGC 10064, gas 
from UGC 10063 could have been stripped off and become 
gravitationally bound to UGC 10064 (which has been seen 
in both simulations and observations, e.g., Barnes & Hern- 
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Figure 15. An optical spectrum of UGC 10063 (SDSS ,1155108.42+254321.3). The inset is a zoom-in of the full spectrum. Narrow 
Balmcr-scries absorption (H/3 through H(J) is marked, and these features yield z = 0.021 ± 0.001, equal to the redshift of UGC 10064 and 
SN 2009dc. Note that Ha is undetected at an expected wavelength of ~6700 A; perhaps Ha emission from residual H II regions fills in 
the absorption line. 



quist 1991; Falcon-Barroso ct al. 2004). This newly acquired 
gas (some of which could be at large galactocentric radii) 
would have likely undergone a burst of star formation soon 
after the two galaxies had their closest encounter (Mihos & 
Hernquist 1994; Barnes & Hernquist 1996), and this might 
be how and when the progenitor of SN 2009dc formed. Both 
SN 2007if and SN 2003fg were also found in hosts with recent 
star formation (Scalzo et al. 2010); in fact, SN 2003fg was 
discovered in a tidal feature (Howell et al. 2006) reminiscent 
of SN 2009dc and UGC 10063. 



4 DISCUSSION 

4.1 Bolometric Luminosity 

To estimate the bolometric luminosity of SN 2009dc, we fol- 
low the method outlined by Howell et al. (2009). Using spec- 
tra of SN 2009dc presented here and spectra of the similar 



SC SN la candidate SN 2007if from Scalzo et al. (2010), 
we take each simultaneous epoch of BVRI photometry and 
warp the spectrum closest in phase to match the photom- 
etry using a smooth third-order spline. We then integrate 
the warped spectrum over wavelengths of 4000-8800 A (i.e., 
from the blue edge of the B-band filter to the red edge of 
/-band filter.) We convert the flux into a luminosity using 
the luminosity distance (di = 97 ± 7 Mpc) adopted in §3.6. 

The above method will only give us the optical bolo- 
metric luminosity. To derive the total bolometric luminosity 
from UV to IR, wc use bolometric corrections from Wang 
et al. (2009b), who found that BVRI only accounts for 
~ 60% of the total luminosity in the case of SN 2005cf, 
a normal SN la. They also found that an additional 20% is 
emitted in the UV and U band, and the final 20% in the 
near-lR. We test to see if the tabulated corrections for the 
UV and U band are reasonable for SN 2009dc by first calcu- 
lating the bolometric luminosity using optical and UV pho- 
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Table 7. Bolometric luminosity for various reddenings 



50 100 150 

Days Past Maximum Light 

Figure 16. Bolometric light curve of SN 2009dc assuming various 
levels of extinction, in comparison with that of SN 2005cf. The 
curve plotted with empty circles assumes no host-galaxy extinc- 
tion, the half-filled circles assume E(B — V)host = 0.1 mag with 
Ry = 3.1, and the filled circles assume E(B — V)host = 0-3 mag 
with Ry =3.1. For comparison, we include the bolometric light 
curve of SN 2005cf as a dashed line (Wang et al. 2009b). We 
indicate our calculated range of values of t-|-i/2 (i-e., the num- 
ber of days after maximum bolometric luminosity when the bolo- 
metric luminosity has decreased to half its maximum value) for 
SN 2009dc within the shaded grey region bounded by dotted lines. 

1/2 = 20.5 days is calculated using a host 
= 0.3 mag and the upper limit of 



The lower bound of tj 



reddening of E{B ~ V)-t, 
t_|_i/2 = 22.8 days assumes no host-galaxy extinction. 



tometry at i = days (using our optical ground-based pho- 
tometry combined with UVOT photometry of SN 2009dc). 
We find that our results arc consistent with a correction 
of 20%. Yamanaka et al. (2009) found that their near-IR 
data for SN 2009dc represented ~20% of the bolometric lu- 
minosity. With this, we cautiously adopt a 40% correction 
for our bolometric luminosity using our well-sampled BVRI 
photometry. Overall, SN 2005cf is a much different SN than 
SN 2009dc, so we include a systematic uncertainty of 20% 
for our bolometric luminosity. 

Figure 16 shows our final bolometric light curves for 
the various levels of host-galaxy extinction we have adopted 
throughout this paper. We also plot the bolometric light 
curve of SN 2005cf for comparison. 

If we use our fiducial nonzero value for the host-galaxy 
reddening, E{B — V')host = 0.1 mag, and our largest reason- 
able value for the host reddening, E{B — V^)host ~ 0.3 mag, 
then we get My = -20.16 ± 0.1 mag (Lboi = (3.3 ± 0.6) x 
10*^ erg s"^) and My = -20.82 ± 0.10 mag (Lboi = (7.3 ± 
1.0) X 10*^ erg s~^), respectively. We note that assuming 
no host-galaxy extinction wo find My — —19.83 ±0.10 mag 
(Lboi = (2.4 ± 0.5) X lO'*^ erg s~^). A summary of these 
results can be found in Table 7. Our range of values of Lboi 
for SN 2009dc matches quite well that of Yamanaka et al. 



E{B - V)host 


-^bol.max 




(mag) 


(10*3 erg s-i) 


(mag) 


0.0 


2.4 ±0.5 


-19.83 ±0.10 


0.1 


3.3 ±0.6 


-20.16 ±0.10 


0.3 


7.3 ± 1.0 


-20.82 ±0.10 



(2009): (2.1-3.3) x 10^^ erg s"^ The most significant differ- 
ence between the two analyses comes from using a different 
range of values ioi E{B — V)host and Ry. 



4.2 ^''Ni Masses and Energetics 

4.2.1 ^^Ni Masses from Amett 's Law 

The vast majority of the bolometric luminosity of SNe la 
comes from the decay of ''^Ni to ^®Co (and subsequently the 
decay of '^'^Co to '"Fe), so we can calculate how much ®®Ni 
was created in SN 2009dc by using "Arnett's law", which 
asserts that the luminosity of a SN la at maximum light is 
proportional to the instantaneous rate of radioactive decay 
of '^^Ni (e.g., Arnett 1982; Stritzinger & Leibundgut 2005). 
Arnett's law is often written 

Lhoi 



Mm = 



aS{tR) 



(1) 



where Mm is the mass of *®Ni present in the ejecta, Lboi is 
the bolometric luminosity at maximum light, a is the ratio 
of bolometric to radioactive luminosities (which is of order 
unity), and S{tR) is the radioactive luminosity per mass of 
®®Ni as a function of the rise time, t_R. 

For SN 2009dc, we have tight constraints on the rise 
time from our photometric data and, as described in §3.1, 
we find a value of 23 ± 2 days. Substituting this into the 
equation for S(tu) from Stritzinger & Leibundgut (2005), we 
get S{tR) « (1.65 ± 0.13) X lO'*^ erg s'^ Mg-^ Stritzinger 
& Leibundgut (2005) find a = 1.2 ± 0.2, which we use for 
our analysis, though we note that this is an upper limit 
(which leads to a lower limit on the ''^'Ni mass) since any 
^®Ni that is above the photosphere will not contribute to the 
SN luminosity. 

Using these values, along with our most likely estimate 
of the bolometric luminosity mentioned above, we calculate 
that SN 2009dc synthesised about 1.7 ± 0.4 Mq of '^'^Ni. If 
we adopt the bolometric luminosity assuming no host-galaxy 
reddening and our largest reasonable host-galaxy reddening, 
along with our calculated rise time and a = 1.2, we get 
~1.2 Mq and ~3.7 M® of ^^Ni, respectively. 



4.2.2 ^^Ni Masses from Late-Time Photometry 

The detection of SN 2009dc in late-time images at levels 
similar to those of normal SNe la offers a stark contrast to 
SN 2006gz, another SC SN la candidate. Late-time pho- 
tometry obtained by Maeda et al. (2009) ~ 360 days af- 
ter maximum found that SN 2006gz had faded significantly 
faster than normal SNe la, casting doubt on the amount 
of ^®Ni synthesised in the explosion. The authors constrain 
Mm < 0.3 Mq assuming that the positrons are nearly com- 
pletely trapped in the ejecta. This is at odds with the ~1 M© 
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of ^^Ni calculated to power the light curve at early times 
(Hicken et al. 2007).'^ To reconcile this difference in nickel 
mass, it was suggested that perhaps SN 2006gz was a highly 
asymmetric explosion, and due to the exact viewing angle 
more ■''^Ni was inferred at early times than was actually 
present (Maeda et al. 2009). However, SN 2009dc appears to 
be spherically symmetric at early times (Tanaka et al. 2010). 
Furthermore, Maeda et al. (2009) posit that dust might have 
formed in the ejecta of SN 2006gz by this time, which would 
reprocess much of the optical light to longer wavelengths. 
However, we do not see strong evidence for dust formation 
in SN 2009dc since the spectral features all appear to evolve 
symmetrically. 

To constrain the nickel mass produced by SN 2009dc 
from late-time photometry, we write out the SN luminos- 
ity as a function of time (e.g., Clocchiatti & Wheeler 1997; 
Maeda et al. 2003, 2009), 

m = Mm (e,(l - e-^(*') + ) e-*/i"-3, (2) 

where t is the number of days from the explosion date, 
Mni is the initial amount of ®®Ni synthesised, £-y = 6.8 x 
10® erg s^^ g^^, T{t) is the time-dependent optical depth to 
7-rays, is the fraction of positrons trapped in the ejecta, 

= 2.4 X 10*^ erg s~^ g^^, and 111.3 days is the e-folding 
time of ^®Co decay. Evaluating r (t) requires defining a model 
with estimates for the kinetic energy of the explosion and 
the total mass of the white dwarf. 

We can place constraints on the ^^Ni mass by only con- 
sidering the luminosity powered by the trapped positrons; 
to do this we set r(t) = in Eqn. (2). Using only the opti- 
cal contribution to the bolometric luminosity over the range 
4000-8000 A from our Keck data taken 281 days past max- 
imum (which is 304 days after explosion), we calculate the 
^^Ni mass for a number of reddening and /g+ values which 
can be found in Table 8. Our calculation of the bolometric 
luminosity does not include a bolometric correction for any 
possible contribution to the bolometric luminosity in the 
NIR. For our nominal values of E{B — F)host = 0.1 mag, 
Rv = 3.1, and /e+ = 1, we calculate a ®®Ni mass of 
1.4 ± 0.3 Mq, consistent with the ^^Ni mass derived from 
Arnett's law and our data at Bmax (§4.2.1). 

Wc also (ictcrmine the ''^'Ni mass using our epoch of 
photometry from 403 days past ma:ximum (which is 426 days 
after explosion). However, without a spectrum to model the 
spectral energy distribution at this phase, wc had to use 
our last spectrum taken 281 days past Bmax, which could 
lead to an inaccurate measure of the bolometric luminosity. 
Furthermore, the detection of the SN was marginal in all 
bands. Using our nominal values of E{B — X^)host =0.1 
mag, Rv = 3.1, and = 1, we calculate a ^®Ni mass of 
0.34 ± 0.06 Mq, clearly much lower than previous estimates 
of the "'^'Ni mass. Wc summarise ^®Ni masses for a range of 
parameters in Table 8. 

The drop in derived ^^Ni mass and expected luminos- 
ity could indicate the formation of dust between t = 281 
and 426 days past maximum light, although this cannot be 



^ Note that this value was derived by assuming a high value for 
the host-galaxy reddening of SN 2006gz, which might have been 
an overestimate. 



Table 8. Ni mass estimates from late-time photometry 



Days past 




E(B - V)hr.„t 


bol 


56 Ni 


Explo s io n 






(■10"° ere s''^) 




304 


0.75 





3.1 ±0.7 


1.4 ±0.3 


304 


1 





3.1 ±0.7 


1.0 ±0.2 


304 


0.75 


0.1 


4.2 ± 0.8 


1.8 ±0.3 


304 


1 


0.1 


4.2 ± 0.8 


1.4 ± 0.3 


304 


0.75 


0.3 


7.9 ± 1.5 


3.4 ± 0.6 


304 


1 


0.3 


7.9 ± 1.5 


2.5 ±0.6 


426 


0.75 





0.26 ±0.05 


0.33 ± 0.06 


426 


1 





0.26 ±0.05 


0.25 ±0.05 


426 


0.75 


0.1 


0.35 ± 0.07 


0.45 ± 0.08 


426 


1 


0.1 


0.35 ±0.07 


0.34 ± 0.06 


426 


0.75 


0.3 


0.64 ±0.12 


0.82 ± 0.15 


426 


1 


0.3 


0.64 ±0.12 


0.62 ±0.12 



^Bolometric luminosity only contains contribution over the range 
4000-8800 A. 



substantiated without IR data or a relatively high-S/N spec- 
trum. The unexpected drop could also be an indication that 

the positron trapping fraction at this phase is significantly 
lower than the nominal value of /e+ = 1. To produce the 
~1.4 Mq of ®®Ni found using our photometry from 281 days 
past maximum, we would need « 0.25. SN 2009dc could 
have undergone an "IR catastrophe" which would cause 
most of the emission to escape at IR wavelengths (as op- 
posed to the optical; Axelrod 1988) as the temperature of 
the SN cools to ~1000 K. This would require using a bolo- 
metric correction to our derived optical bolometric luminos- 
ity. Alternatively, the low derived ^''Ni mass could suggest 
that the late-time light curve of SN 2009dc is not powered 
by the decay of ^®Co. This conclusion, however, is at odds 
with the ^®Ni mass derived from our epoch of photometry 
from t — 281 days past maximum light. Without ncar-IR 
data or a spectrum at this epoch, it is difficult to reconcile 
these differences. We again caution that the spectral energy 
distribution is ill constrained at this phase. 



4-2.3 WD Masses and Kinetic Energies 

Using the technique developed by Howell et al. (2006) , we at- 
tempt to combine our observations of SN 2009dc with prop- 
erties of WDs and SNe la in order to calculate the mass of 
the progenitor WD as well as the amount of ^^Ni produced 
in the explosion. We note that this method is independent 
of the one used in the previous section which utilitized the 
late-time photometry of SN 2009dc. As will be mentioned 
below, there are a few assumptions and simplifications that 
go into this analysis, and thus the final values should prob- 
ably not be taken literally. However, the sense of the rela- 
tionships and trends these calculations indicate should be 
reliable and will show that SN 2009dc likely had a SC WD 
progenitor. 

The ejecta of SNc la derive their kinetic energy (Ek) 
from the energy released during explosive nucleosynthesis 
(En), but the ejecta must first overcome the binding energy 
(Eb) of the WD progenitor (Branch 1992). Therefore, we can 
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write 

Ek = — — = En — Ed, (3) 

where Mwd is the mass of the WD just before explosion and 
w is a representative velocity of the ejecta. A single value for 
velocity only relates to the kinetic energy of the SN averaged 
over the entire ejecta, and thus there has been disagreement 
as to which actual value should be used for a given SN. 
Previous studies have used the velocity of the Si II A6355 
feature at 10 and 40 days past maximum (Benetti et al. 
2005; Howell et al. 2006). 

Since about 70% of IGEs produced in SNe la are in the 
form of ^®Ni (Nomoto et al. 1984; Khokhlov et al. 1993), the 
total mass of the WD can be related to the nickel mass by 

Mm = 0.7Mwd/ige, (4) 

where /ige is the fractional composition of IGEs in the SN 
ejecta. However, the exact amount of IGEs that are made 
up of ^®Ni may not be the same for all SNe la and likely 
dependent on other, external factors such as metallicity (e.g., 
Howell et al. 2009, and references therein). 

Furthermore, Branch (1992) found that fusing equal 
parts carbon and oxygen all the way to iron yields epe = 
1.55 X 10^^ erg M©^^, and burning the same material only 
up to silicon liberates about 76% of that energy. Thus we 
can write 

En = EFeMwD (/iGE + 0.76/ime) , (5) 

where /ime is the fractional composition of IMEs in the 
ejecta. 

Finally, we define the fractional composition of un- 
burned carbon in the SN ejecta as /c.* Based on our defi- 
nitions we can write 

1 = /c + /ime + /ige. (6) 

Putting all of this together, we relate various properties of 
the progenitor WD to properties of the SN ejecta it created 
using 

However, given the assumptions and simplifications that 
went into this "equation" , wo stress that it should be used 
more as an illustration of how the different parameters re- 
late to each other and less as a tool to actually calculate 
WD and ^''Ni masses. 

Values of Ei, range from about 0.5 x 10 '^ erg for a WD 
with a mass of 1.4 Mq to about 1.3 x 10^^ erg for a WD 
with a mass of 2 M© and a central density of 4 x 10^ g cm~^ 
(Yoon & Langer 2005). We perform separate calculations 
for each pair of these values of Eb and Mwd- Tanaka et al. 
(2010) measure the Si II A6355 feature to be blueshifted 
by 7200 km s~^ in their spectrum of SN 2009dc obtained 
6 days past maximum, while we find the same feature to be 
blueshifted by 5000-6000 km s"'^ in our spectra taken 35, 
52, and 64 days past maximum. Since v has been evaluated 

® This parameter does not appear In the equation for E„ ; it rep- 
resents the amount of unbumed material, and thus does not con- 
tribute to the energy released by nucleosynthesis. 



at 10 and 40 days past maximum for similar analyses pre- 
viously (Benetti et al. 2005; Howell et al. 2006), we will use 
a range of velocities from 7000 km s~^ to 5500 km s~^ for 
our analysis. Finally, since there is clear indication of un- 
burned material in the ejecta of SN 2009dc (see §3.3.2), we 
know that fc must be nonzero. However, based on models 
of SN la ejecta (e.g., Nomoto et al. 1984) and the relative 
scarcity of carbon detections in SNe la (Marion et al. 2006; 
Thomas et al. 2007), we will only use values in the range 
0.05-0.3 for fc. 

With such large uncertainties and so many parameters, 
we will clearly calculate a huge range of ®®Ni masses, and we 
realize that a few assumptions and model-dependent values 
have come into the derivation. However, this is still a useful 
line of reasoning in our attempt to determine both the mass 
of the WD progenitor and the amount of ^^Ni synthesised 
by SN 2009dc. We therefore forge ahead and attempt to 
eliminate at least some of the most extreme situations. 

Almost all of parameter space is ruled out when we use 
a 1.4 M0 WD (and its associated binding energy), since we 
calculate negative nickel masses for these cases. The largest 
^^Ni mass we can reasonably derive for a 1.4 Mq WD (using 
fc = 0.3 and v = 7000 km s^^) is 0.06 M0, nearly an order 
of magnitude smaller than the amount of nickel produced by 
a normal SN la (e.g., Nomoto et al. 1984; Kasen et al. 2008). 
We have also shown that SN 2009dc is much more luminous 
than the average SN la. Therefore, if it did have a 1.4 
progenitor, then we would need to invoke a large, hitherto 
unknown energy source to power its light curve, which seems 
highly unlikely. 

The situation changes when we use a 2 M© WD (with 
a central density of 4 x lO'"* g cm""^), put forth as a possible 
SC SN la progenitor by Yoon & Langer (2005). Again, a 
significant part of parameter space is eliminated, since at 
low values of fc we calculate negative (or extremely small) 
nickel masses. Assuming that /ime must be > 0.1 (given that 
we detect IMEs in the spectra for months after maximum; 
see, e.g., Nomoto et al. 1984), we calculate nickel masses 
as high as ~ 1 Mq. In order to match our nominal ®^Ni 
mass range of 1.4-1.7 M© and satisfy our constraints on the 
fractional composition of elements mentioned previously, we 
must resort to using a WD progenitor more massive than 
2M0. 

Again, the final values calculated above are almost cer- 
tainly not the true answer. However, the range of plausible 
values ijresented, as well as the relationships between the 
parameters of the SN and its WD progenitor, yet again in- 
dicate that the progenitor of SN 2009dc was Ukely a SC WD. 

The various analyses presented above seem to strongly 
favor the conclusion that the progenitor of SN 2009dc was a 
SC WD with a mass of probably greater than ~2 Mq . These 
analyses also indicate that SN 2009dc most likely produced 
1.4-1.7 Mq of ^®Ni (assuming our nominal value for the 
host-galaxy reddening and our peak bolometric luminosity). 
More than 1 Mq of ^^Ni was almost certainly created by 
SN 2009dc, and the actual value couldhe as high as ~3.3 Mq. 
This matches well with the conclusions of Yamanaka et al. 
(2009), who calculate a similar range of nickel masses for 
SN 2009dc (1.2-1.8 Mq). Part of the difference can be ac- 
counted for by the fact that the two studies use different 
values of host-galaxy reddening (causing the derived bolo- 
metric luminosities to differ somewhat; see §4.1). However, 
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the majority of the difTerence in derived ^^Ni masses comes 
from the assumed rise time of SN 2009dc. Yamanaka et al. 
(2009) adopt a rise time of 20 days based on comparisons to 
typical SNe la and SN 2006gz, while we use a rise time of 
23 days based on our pro-maximum photometry. The longer 
rise time used in our study leads to a larger derived nickel 
mass for SN 2009dc. 

SN 2003fg had a derived progenitor mass of ~2.1 M© 
and produced ~1.3 M© of ^®Ni (Howell et al. 2006), and 
these values are quite similar to those we calculate for 
SN 2009dc. SN 2006gz was estimated to have produced 1- 

1.2 Mq of ^^Ni, which is on the low end of our range for 
SN 2009dc, and it was also claimed to have a SC WD pro- 
genitor (though no attempt was made to further constrain 
the progenitor mass; Hicken et al. 2007). 

4.3 Comparison to Theoretical Models 

Any theoretical model which is postulated to explain 
SN 2009dc, with or without a SC WD, must be able to re- 
produce the observed peculiarities for which we have very 
tight constraints: (a) high luminosity even when assuming 
no host-galaxy reddening, (b) relatively long light-curve rise 
time, (c) relatively slow photometric decline and late-time 
spectroscopic evolution, (d) the presence of carbon in spec- 
tra near maximum brightness, (e) the presence of silicon in 
spectra as late as a few months past maximum brightness, {fj 
IGEs dominating the spectra at late times, and (g) mostly 
spherically symmetric ejecta near maximum with possible 
clumpy layers of IMEs (Tanaka et al. 2010). 

Below we consider both SC and non-SC models, all of 
which involve the thermonuclear explosion of a WD. How- 
ever, the similarities between SN 2009dc and SN 2002cx and 
its brethren (see §3.4.5) may argue that the progenitor and 
explosion scenario for these objects are all linked. Could 
SN 2009dc be a SN 2002cx-like object but with a much 
higher luminosity? Perhaps, although attempts to explain 
the origins of even the best-observed SN 2002cx-like objects 
are still ongoing (e.g., Jha et al. 2006; Phillips et al. 2007; 
Valenti et al. 2009; Li et al. 2010). 

4.3.1 SC Models 

Two- and three-dimensional models of differentially rotat- 
ing massive WDs have been presented in the literature, 
and calculations show that SC WDs with masses as large 
as 2 Mq are possible by accretion from a nondegenerate 
companion, while masses up to ~1.5 M© are possible from 
a double-degenerate merger (Yoon & Langer 2005; Yoon 
et al. 2007, respectively). More recent studies of double- 
degenerate mergers have even shown possible SN la progen- 
itors with total masses approaching 2.4 Mq (Greggio 2010). 
This is encouraging for the case of SN 2009dc, since our en- 
ergetics arguments imply that its progenitor is likely greater 
than ~2 Mq. However, calculations by Piro (2008) suggest 
that differential rotation is unlikely for WDs accreting from 
a nondegenerate companion, and thus their masses cannot 
exceed the Chandrasekhar mass by more than a few percent. 

Yoon & Langer (2005) also consider how much '^Ni 
would be produced by a SN la whose progenitor is ~2 Mq, 
calculating values of 0.4-1.3 Mq. This large range does en- 
compass the lower end of our range of ®®Ni for SN 2009dc. 



However, only specific kinetic energies (i.e., kinetic ener- 
gies per unit mass) which are lower than what we find for 
SN 2009dc by a factor of ~2 were used (Yoon &: Langer 
2005). Adopting lower specific kinetic energies will tend to 
decrease the nickel yield for a given WD mass, according to 
their models. It is possible that increasing the specific ki- 
netic energies used in the models of Yoon & Langer (2005) 
to values that better match SN 2009dc will in fact reproduce 
our derived nickel mass, but this may require pushing the 
models into regimes where they are no longer valid. 

Pfannes et al. (2010) use aspects of the WD models of 
Yoon & Langer (2005) and simulate prompt detonations of 
massive, rapidly rotating WDs. None of their models ap- 
pears to reproduce the fact that we detect unburned fuel 
(i.e., carbon) and IMEs near maximum brightness at low 
velocities, and subsequently IMEs and IGEs at later times 
(also at low velocities). Even more at odds with SN 2009dc, 
Pfannes et al. (2010) predict that the IMEs and unburned 
material should have similar spatial distributions in the 
ejecta (specifically, in a torus in the equatorial plane), but 
this seems unlikely given the strong line polarizations seen 
by Tanaka et al. (2010) in Si II and Ca II features and the 
lack thereof in two different C II features. Thus, we con- 
clude that the models of Pfannes et al. (2010) do not seem 
to reproduce the observations of SN 2009dc. 

Models of SC WDs and their evolution in close binary 
systems with nondegenerate companions can also be found 
in the literature (Chen &: Li 2009). These evolutionary cal- 
culations include the effects of varying the orbital period of 
the binary systems, the metallicity and mass-transfer rate 
of the binary companion, and (most importantly) the mass 
of the WD. Chen & Li (2009) find that WDs with initial 
masses of ^1.2 M,:,,, midcr the right mass-transfer condi- 
tions, can accrete up to masses of about 1.4-1.8 M© before 
exploding as a SN la (though they find that most of these 
WDs explode with masses not much above 1.4 Mq). Thus it 
seems somewhat unlikely that these models can explain the 
progenitor we suggest for SN 2009dc; even if the occasional 
WD can accrete up to 1.8 Mq, this is at the lowest end of 
our range of WD progenitor masses. 

Spherically symmetric, one-dimensional radiation 
transport calculations for normal and SC WDs have been 
carried out by Maeda & Iwamoto (2009). They find that 
based on the light-curve shapes, photospheric velocities, 
peak bolometric luminosities, and peak effective temper- 
atures, SN 2006gz likely came from a SC WD'' while 
SN 2003fg did not. Since the observables of SN 2003fg 
used for their analysis have very similar values to those 
of SN 2009dc, it would seem that their models imply that 
SN 2009dc is also not a SC SN la. Maeda & Iwamoto (2009) 
point out that SN 2003fg (and thus probably SN 2009dc as 
well) could be a SC SN la if the progenitor star were highly 
aspherical, but again this seems unlikely for SN 2009dc 
from the spectropolarimetric data (Tanaka et al. 2010). 

The primary independent variable used in the calcu- 
lations of Maeda & Iwamoto (2009) is the number 
of days after ma^cimum bolometric luminosity when the 



^ Maeda & Iwamoto (2009) claim that most of the emission likely 
shifted to the near-IR and mid-IR at late times in order to explain 
the relatively faint late-time observations of Maeda et al. (2009). 
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bolometric luminosity has decreased to half its maximum 
value (Contardo et al. 2000). For SN 2006gz they measure 
t+i/2 = 18 days from the publicly available photometry, 
but for SN 2003fg they convert the stretch value published 
by HowoU et al. (2006) to Amu and then to a t+i/2 value 
of 13.5 days (adopting an empirical linear fit to Amis and 
t+i/2 values from Contardo et al. 2000). Using their conver- 
sions between Amis and t+i/2, we calculate t+i/2 14 days 
(which, unsurprisingly, is nearly the same as the value calcu- 
lated for SN 2003fg by Macda & Iwamoto 2009). However, 
when we measure t+i/2 directly from our light curve, we 
get a minimum value of about 20.5 days (when we use our 
maximum plausible host-galaxy reddening), which is nearly 
50% larger! Furthermore, if one converts the Amis value 
of SN 2006gz to t+i/2 using the linear fit, one again finds 
t+i/2 ~ 14 days, suggesting that the linear conversion de- 
rived by Maeda & Iwamoto (2009) is not accurate for such 
low values of Amis. This is not wholly unexpected since the 
Amis values of SNe 2006gz, 2003fg, and 2009dc (or t+i/2 
as measured from light curves) are all below (or above) the 
values used to derive the linear fit (Contardo et al. 2000). 
In addition, the conversion from stretch to Amis used by 
Maeda & Iwamoto (2009) for SN 2003fg did not take into 
account the fact that the definition of stretch has evolved 
as new SN la light-curve templates have been constructed 
(e.g., Conley et al. 2008). 

We can now compare SN 2009dc to the analysis of 
Maeda & Iwamoto (2009) using the actual values of t+i/2 ~ 
20.5-22.8 days. Given the large bolometric luminosity and 
t+i/2 value, as well as the extremely low photospheric veloc- 
ity of SN 2009dc, none of their models (using normal or SC 
WDs) appears to be viable. Their "normal WD" model with 
MwD = 1.39 M0 and Mm =0.6 M3 accounts for the low 
velocity with large t+i/2, but it underpredicts Lboi by a fac- 
tor of a few. Some of the SC WD models (Mwd ^ 2.0 Mq 
and Mni = 1.0 M©) of Maeda & Iwamoto (2009) can al- 
most account for the large values of Lboi and t+i/2 seen in 
SN 2009dc, but they then ovcrpredict the photospheric ve- 
locity of SN 2009dc by a factor of 2 or so (with the models 
having the highest mass WDs coming the closest to match- 
ing the observed velocities). 

While no single model of Maeda & Iwamoto (2009) 
clearly reproduces our observations of SN 2009dc, its prop- 
erties seem to be on the outskirts of the parameter space 
explored by their analysis. This hints at the possibility that 
more extreme WD and ''^Ni masses may be required to 
match the observational data of SN 2009dc. In addition, it 
is possible that multi- dimensional analyses are required to 
truly capture the underlying physics of a SC WD explosion. 

Furthermore, a few synthetic late-time spectra of 
SC SNe la have been presented in the literature (e.g., 
Maeda et al. 2009). Their models take as inputs a WD 
progenitor meiss, a WD central density, and mass frac- 
tions of burning products, and then use a one-dimensional 
Monte Carlo radiation transport code, along with ioniza- 
tion/recombination equilibrium and rate equations, to cal- 
culate synthetic light curves and spectra. Both of their 
SC SN la models roughly match our SN 2009dc photom- 
etry. In addition, their late-time spectra are similar to our 
day 281 spectrum of SN 2009dc. Maeda et al. (2009) note 
that as they increase their models' progenitor mass, the 
[Fe II] /[Ca II] feature near 7200 A becomes stronger rela- 



tive to the blends of [Fe II] and [Fe III] near 3800-5500 A. 
This is seen in SN 2009dc (as compared to SNe la with more 
normal peak luminosity) , and it further supports our finding 
that the progenitor of SN 2009dc was a SC WD. However, 
it should be noted that the SC SN la models of Maeda et al. 

(2009) only use WDs with masses of 2 and 3 Mq, and ^®Ni 
masses of only 1 M©, which is on the low end of our range 
of calculated values for the ^®Ni yield of SN 2009dc. 

4.3.2 Non-SC Models 

A number of models which employ Chandrasekhar-mass 
WDs as the progenitors of super- luminous SNe la (ones that 
we would consider possibly SC SNe la) have been proposed 
(HiUebrandt et al. 2007; Sim et al. 2007; Kasen et al. 2008, 
2009). These models usually invoke the off-centre ignition 
of a normal WD progenitor which leads to nuclear burn- 
ing (and thus '''Ni production) that is peaked away from 
the centre of the WD. This off-centre nickel blob would in- 
crease the observed luminosity if the blob were offset from 
the centre of the WD toward the observer, with the maxi- 
mum effect occurring when it is offset directly along the line 
of sight. These viewing angles also lead to the fastest light- 
curve rise times in such simulations (Sim et al. 2007; Kasen 
et al. 2008). 

Once again, we can compare our observed values for 
SN 2009dc to the models. Since the nickel blob is offset from 
the centre, asphericity is introduced into the explosion by 
construction (Kasen et al. 2008). If there is in fact a blob 
of ®®Ni in SN 2009dc and it is offset from the progenitor's 
centre directly along our line of sight (thus maximizing the 
measured luminosity), then perhaps there is still azimuthal 
symmetry in the explosion. This may account for the low 
levels of continuum polarization measured by Tanaka et al. 

(2010) , as well as the higher levels of IME line polarization, 
but it seems tenuous at best. 

Observing along the axis of the offset blob of ^^Ni, HiUe- 
brandt et al. (2007) produce light curves that pealc at a bolo- 
metric magnitude of Mboi ^ —19.9 mag, Kasen et al. (2008) 
show light curves that get as bright as Mb ~ —20 mag, and 
Kasen et al. (2009) claim models with luminosities as high 
as 2.1 X lO**"^ erg s~^. These brightest magnitudes and lu- 
minosities that can be obtained by off-centre explosions are 
nearly equal to our lower limits for SN 2009dc (assuming no 
host-galaxy reddening), and are likely below the true values. 
Furthermore, we note that the maximum amount of ^®Ni 
obtained from these explosions is 0.9-1.1 M© (HiUebrandt 
et al. 2007; Kasen et al. 2009), which is once again at the 
lowest end of our range of calculated values for SN 2009dc. 

Finally, as mentioned above, the viewing angles that 
maximize the observed peak magnitudes also minvrmze the 
rise time. We find a relatively long rise time of 23 ± 2 days 
for SN 2009dc, and our first detection of the SN implies that 
its rise time must be >21 days. This is significantly longer 
than the rise times for these viewing angles as derived from 
the models (~12 days and ~18 days; HiUebrandt et al. 2007; 
Kasen et al. 2008, respectively). Thus, it seems that none of 
these models which include a Chandrasekhar-mass WD is 
viable for SN 2009dc. 

It should be mentioned, however, that these models 
mostly assume expansion velocities of "normal" SNe la (e.g. 
Sim et al. 2007) and we have shown that the expansion veloc- 
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ity of SN 2009dc is significantly lower. This is evident in the 
synthetic spectra derived from these models; from all view- 
ing angles, they resemble early-time normal SN la spectra 
much more than the spectra of SN 2009dc near maximum 
brightness (Kascn ct al. 2008). Thus, perhaps it is not too 
surprising that these specific examples of non-SC models do 
not match the observations of SN 2009dc. 



5 CONCLUSIONS 

In this paper we have presented and analysed optical pho- 
tometry and spectra of SN 2009dc and SN 2007if, both of 
which are possibly SC SNe la. Our photometric and spectral 
data on SN 2009dc constitute one of the richest datasets ever 
published on a SC SN la candidate. Our well-sampled light 
curve follows SN 2009dc from about 1 week before maximum 
brightness until about 5 months past maximum, and shows 
that SN 2009dc is one of the slowest photometrically evolv- 
ing SNe la ever observed. We derive a rise time of 23 days 
and Ami5(B) = 0.72 mag, which are two of the most ex- 
treme values for these parameters ever seen in a SN la. As- 
suming no host-galaxy reddening, we derive a peak bolo- 
metric luminosity of about 2.4 x 10''^ erg s~^, though this is 
almost certainly an underestimate since we observe strong 
evidence for at lease some host reddening. Using our nonzero 
values for E{B — V^)host, the peak bolometric luminosity in- 
creases by about 40%-200%. 

Spectroscopically, SN 2009dc also evolves relatively 
slowly. Strong C II absorption features (which are rarely 
observed in SNe la) arc soon in the spectra near maximum 
brightness, implying a significant amount of unburnod fuel 
from the progenitor WD in the outer layers of the SN ojccta. 
Si II absorption also appears in our spectra of SN 2009dc 
and remains visible even 2 months past maximum. Our 
post-maximum spectra arc dominated by a forest of IGE 
features and, interestingly, resemble spectra of the peculiar 
SN la 2002cx. Finally, the spectra of SN 2009dc all show very 
low expansion velocities at all layers (i.e., unburned carbon, 
IMEs, and IGEs) as compared to other SNe la. This may be 
explained by a massive WD progenitor which consequently 
has a large binding energy. Even though the expansion ve- 
locities are small, we see no strong evidence in SN 2009dc for 
a velocity "plateau" near maximum light like the one seen 
in SN 2007if (Scalzo et al. 2010). 

Using various luminosity and energy arguments, we cal- 
culate that the progenitor of SN 2009dc is possibly a SC WD 
with a mass greater than ~2 M©, and that at least ~1 Mq 
of ''''Ni was likely formed in the explosion (though the most 
probable value is in the range 1.4-1.7 Mq). These values 
are larger than (or about as large as) those calculated for 
any other SN la ever observed. We propose that the host 
galaxy of SN 2009dc underwent a gravitational interaction 
with a nearby galaxy (UGC 10063) in the relatively recent 
past, and that this could have induced a sudden burst of 
star formation which may have given rise to the progeni- 
tor of SN 2009dc and turned UGC 10063 into the "post- 
starburst" galaxy that we observe today. We also compare 
our observed quantities for SN 2009dc to theoretical models, 
and while no model seems to match or explain every aspect 
of SN 2009dc, simulations show that SC WDs with masses 
near what we calculate for the progenitor of SN 2009dc can 



possibly form, likely from the merger of two WDs. Further- 
more, models of extremely luminous SNe la which employ 
a Chandrasekhar-mass WD progenitor cannot explain our 
observations of SN 2009dc. 

Thus, taking all of these extreme values into account, 
we conclude that SN 2009dc is very likely a SC SN la. As 
mentioned previously, many of the observed peculiarities of 
SN 2009dc are also seen in SN 2003fg and SN 2007if. There- 
fore, we concur with Howell et al. (2006) and Scalzo et al. 
(2010) that both SN 2003fg and SN 2007if (respectively) are 
also probably SC SNe la. However, given their fairly normal 
expansion velocities and relative weakness (or even absence) 
of C II features near maximum brightness, it seems that 
SN 2006gz and SN 2004gu are less likely to be SC SNe la. 

New large transient searches such as Pan-STARRS 
(Kaiser et al. 2002) and the Palomar Transient Factory (Rau 
et al. 2009; Law et al. 2009) will probably find many SC or 
other super-luminous SNe la in the near future. Since it 
seems that they cannot be standardized in the same way as 
most SNe la, they will need to be handled separately or ig- 
nored in cosmological surveys which will use large numbers 
of SNe la. However, the simulations of Chen & Li (2009) 
show that donor stars with lower metallicities (e.g.. Popula- 
tion II stars) are less likely to form WDs with masses greater 
than 1.7 Mq than higher metallicity stars. Thus, it is pos- 
sible that contamination levels from SC SNe la, which are 
already rare at low redshifts (i.e., average metallicity), may 
be relatively small in medium or high-redshift surveys. 
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